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Abstract 



We present realistic models of pion-exchange three-nucleon interactions 
obtained by fitting the energies of all the 17 bound or narrow states of 
3 < A < 8 nucleons, calculated with less than 2% error using the Green's 
function Monte Carlo method. The models contain two-pion-exchange terms 
due to ttN scattering in S- and P-waves, three-pion-exchange terms due to 
ring diagrams with one A in the intermediate states, and a phenomenological 
repulsive term to take into account relativistic effects, the suppression of the 
two-pion-exchange two-nucleon interaction by the third nucleon, and other ef- 
fects. The models have five parameters, consisting of the strength of the four 
interactions and the short-range cutoff. The 17 fitted energies are insufficient 
to determine all of them uniquely. We consider five models, each having three 
adjustable parameters and assumed values for the other two. They reproduce 
the observed energies with an rms error < 1% when used together with the 
Argonne vi$ two-nucleon interaction. In one of the models the ttN S-wave 
scattering interaction is set to zero; in all others it is assumed to have the 
strength suggested by chiral effective field theory. One of the models also as- 
sumes that the ttN P-wave scattering interaction has the strength suggested 
by effective field theories, and the cutoff is adjusted to fit the data. In all 
other models the cutoff is taken to be the same as in the v is interaction. The 
effect of relativistic boost correction to the two-nucleon interaction on the 
strength of the repulsive three-nucleon interaction is estimated. Many calcu- 
lated properties of A < 8 nuclei, including radii, magnetic dipole and electric 
quadrupole moments, isobaric analog energy differences, etc., are tabulated. 
Results obtained with only Argonne v' 8 and vi$ interactions are also reported. 
In addition, we present results for 7- and 8-body neutron drops in external 
potential wells. 

PACS numbers: 21.10.-k, 21.45.+V, 21.30.+y, 13.75.Cs, 12.40.Qq, 
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I. INTRODUCTION 



One of the primary goals of nuclear physics is to understand the stability, structure, 
and reactions of nuclei as a consequence of the interactions between individual nucleons. 
However, these interactions are not known from first principles; they are modeled with pa- 
rameters to be determined from data. Significant advances have been made during the last 
decade in the ah initio calculation of nuclear properties starting from these realistic models 
of the nuclear force, which allow us to test the predictions of such models with unprece- 
dented accuracy, and to refine them. With our collaborators, we have carried out a series of 
many-body calculations of light nuclei |I|fJ] and nuclear and neutron star matter || using a 
Hamiltonian that contains both two- and three-nucleon potentials. The light nuclei calcula- 
tions use the Green's function Monte Carlo (GFMC) method and have been demonstrated to 
give nuclear binding energies for up to eight-body nuclei with a precision of better than 2%. 
The matter calculations are less accurate but provide important constraints on the Hamil- 
tonian. These calculations have used the Argonne t>ig (AV18) model [[|] of the two-nucleon 
interaction, t^-, and the Urbana IX (UIX) model || of the three-nucleon interaction, Vyk- 

The results for light nuclei are summarized in Fig. [I], where we compare the calculated 
and experimental binding energies for all the ground or narrow, low-lying, excited states 
of nuclei with up to eight nucleons (neglecting isobaric analog states). In addition to the 
predictions of the AV18/UIX Hamiltonian, we show the results (most newly calculated for 
the present paper) of using just the two-body AVI 8 interaction by itself. We see that AVI 8 
alone predicts some key features of nuclear structure correctly, such as the proper ordering 
of excited states and the rapid saturation of the binding above 4 He. However, with the 
exception of 2 H, it underbinds all nuclei, and this failure grows rapidly with increasing A. 
With just the two-nucleon force acting, the Borromean nuclei 6 He and 8 He are not stable 
and the lithium nuclei are only marginally so. 

The addition of the UIX model of Vijk fixes the binding energy of 3 H and 4 He and 
significantly improves the binding of the p-shell nuclei. However, AV18/UIX still underbinds 
as A increases, and also as N — Z increases. In particular, 8 He is more underbound than 8 Be, 
indicating a problem with the isospin dependence of this interaction model. The relative 
stability of the lithium nuclei is improved, but the Borromean helium nuclei are still unbound. 
Additional calculations of wider, higher-lying, excited states not shown in Fig. [I] indicate 
another problem with the AV18/UIX model: the underprediction of spin-orbit splittings 
among spin-orbit partners such as the | and | states in 5 He. 

In this paper we investigate new models of V^k that largely correct these failings and 
give a very good description of the spectrum of light nuclei. Studies of nuclear and neutron 
star matter with these new models will be reported in a separate paper. 

The theory of strong interactions has not yet progressed enough to permit a first- 
principles determination of the two- and three-nucleon interactions with the accuracy re- 
quired to calculate nuclear binding energies. The interactions must be determined phe- 
nomenologically. Modern, realistic models of Vy are obtained by fitting the ~4300 data 
below 350 MeV in the Nijmegen iVTV-scattering data base || with a \ 2 ~ 1 P er degree of 
freedom. The Nijmegen data base is said to be complete, i.e., the included data determine 
all the relevant phase shifts and mixing parameters. Thus v y fitted to it are well determined 
and generally give very similar predictions of the properties of three- and four-body nuclei, 
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as will be discussed below. 

In contrast it is much more difficult to construct realistic models of V^k by simply fitting 
three-nucleon scattering data, which is dominated by the pairwise forces. The number of 
operators that can contribute to V^k is very large, and until recently, the number of observ- 
ables that could both be observed and accurately calculated was small. Recent advances 
in three-nucleon scattering calculations, based on correlated hyperspherical harmonic |7j 
and Faddeev [|S] methods, and in high-precision Nd scattering experiments, hold significant 
promise for testing models of V^k in this regime. However, the binding energies and excita- 
tion spectra of light nuclei also contain a great deal of information, and are in fact the only 
current means to investigate T — | forces. 

An additional concern is that the V^k obtained by fitting nuclear data may depend 
strongly on the model of Vij used in the Hamiltonian. The Vijk will naturally depend 
upon the chosen to some extent. For example, two equivalent models of Vij, related 
by a unitary transformation, will have different but related V^k associated with them ||. 
However, combinations of and Vijk related by unitary transformations will naturally 
predict the same observables. 

Models of V^k based on the elimination of field variables date back to the work of 
Primakoff and Holstein ||10|| . The first modern meson-exchange model for nuclear Vijk was 
proposed by Fujita and Miyazawa (FM) [|ll[|; it contained only the two-pion-exchange three- 
nucleon interaction V 2n,PW due to scattering of the pion being exchanged between two 
nucleons by a third nucleon via the P-wave A-resonance. This interaction is attractive in 
nuclei and nuclear matter. Later theoretical models, such as Tucson-Melbourne (TM) |l2j 
and Brazil ]T3] included the V 2n,sw due to ttN S-wave scattering and V 27T,PW from all P- 
wave scattering. In the recent Texas model these two-pion-exchange contributions to V^k 



have been predicted using chiral symmetry |14j|. The FM and later models have similar 



forms for V 2n,PW , but the predicted strength of the long-range part of V 27T,PW in the later 
models is almost twice that in FM. 

The main failures of a nuclear Hamiltonian containing only two-nucleon interactions 
include the underbinding of light nuclei, as discussed above, and an overestimate of the 
equilibrium density of nuclear matter. An attractive V 2w addresses the first failure while 
making the second worse ||15|| . The Urbana models of Vijk contain only two terms, the V 2n,PW 
and a phenomenological, repulsive V R . The strengths of the two interactions in the most 
recent Urbana model, UIX, were obtained by reproducing the energy of 3 H via a GFMC 
calculation and the density of nuclear matter by approximate variational calculations [^Q. 
The repulsive term V R in Vijk is essential to prevent nuclear matter from being too dense 
and overbound. 

The expectation value of the NNN potential is much smaller than that of the NN po- 
tential in nuclei. For example, the ratio of contributions of the UIX and AV18 potentials 
in A < 8 nuclei is < 0.1 |1]|§. However, the V^k gives a relatively much larger contribution 
to nuclear binding energies due to the significant cancellation between the positive kinetic 
energy and the negative NN potential. It is this feature that allows us extract information 
on V^k by studying the spectrum of light nuclei. 

In the present work we fit the energies of 17 states of up to eight nucleons, calculated 
by the GFMC method with < 2% error, to construct more realistic models of Vijk- These 
models are for use with the AV18 v^. In addition to the already mentioned V 2lT,PW , V 27T,SW , 
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and V R terms, they contain three-pion-exchange rings with A intermediate states, V 37r ' . 
All the terms are static; their spin-isospin and spatial dependence is taken from theoretical 
models, and their strengths are varied to fit the observed energies. 

The new models are refered to as "Illinois" models; five versions, Illinois-1 to -5 (desig- 
nated IL1 to IL5) are presented in this paper. The Hamiltonians using AV18 and these Vijk 
are refered to as AV18/IL1, etc. For each model, two to three of the available five parameters 
were adjusted to fit the binding energies of the 17 states assuming plausible values for the 
other parameters. The IL1 and IL2 models have short-range cutoffs taken from AV18 , while 
IL3 uses the strength predicted by chiral perturbation theory jTJJ for V^jJ, and adjusts the 
cutoff to fit the energies. IL4 and IL5 are further variations of IL2. The qualities of the fits 
are good, and the extracted strength parameters have plausible values. This suggests that 
strengths of additional terms in Vijk cannot be determined from the data included in the 
present work. It is also possible that additional terms in V^k are weaker than V 2n,sw and 
V 3n > AR , which in turn are weaker than the V 27T ' PW and V R considered in the older Urbana 
models. 

Several relativistic effects are contained in the two- and three-nucleon potentials fitted 
to experimental data. However, the boost correction, 5v(Pij), to the two-body interaction 
is omitted in nonrelativistic Hamiltonians containing u^- fitted to the scattering data in the 
two-nucleon center of mass frame. This many-body effect arises from the motion of the 
center of mass of the ij pair of nucleons in the presence of the other nucleons. Initially 
the boost interaction 8v(Pij) is neglected in the GFMC calculations. It is subsequently 
treated as a first-order perturbation. The contribution of 5v(Pij) to the binding energy of 
light nuclei is nearly proportional to that of V . The final value of the strength of V R can 
be adjusted to reproduce the observed energies when the perturbatively computed 5u(Py) 
contribution is included. 

A brief review of the NN interaction, including relativistic corrections, is given in Sec. 
II. The new Illinois models of Vijk are presented in Sec. III. The GFMC calculations of light 
nuclei are briefly described in Sec. IV. The nuclear energies calculated with AV18 and its 
approximation AV8', as well as those including the new Illinois three-nucleon potentials are 
reported in Sec. V. A number of results obtained with the new Illinois models for the light 
nuclei, including proton and neutron distribution radii, magnetic and quadrupole moments, 
and isobaric analog energy differences, are also given in Sec. V. In addition, we report results 
obtained for drops of 7 and 8 neutrons in an external potential well to provide constraints 
for energy-density functionals of neutron-rich nuclei jrj|. Our conclusions are given in Sec. 
VI. 



II. THE TWO-NUCLEON INTERACTION 

We use the Hamiltonian: 

# = £-|^ 2 + £ % + E vm, (2.1) 

i 1 i<j i<j<k 

containing kinetic, two- and three-nucleon interaction energies. The mass difference between 
the proton and the neutron is taken into account by letting be the mass of proton or 
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neutron according to the isospin of nucleon i, and both strong and electromagnetic isovector 



and isotensor terms are included in the i>y. 



The Argonne Vi$ two-nucleon potential [ffl contains v n , the one-pion-exchange potential 
with a short-range cutoff, v R representing all other strong interaction terms, and v 7 , a very 
complete treatment of the electromagnetic interaction: 



% =vL + ^ + v7 i . (2.2) 



It can be expressed as a sum: 



Vij = v p( r ij)°ij > (2-3) 
p 



in which Ofj are operators, and v p (rij) depend only on the interparticle distance r^. The first 
six operators are the only possible isospin-conserving static ones, i.e., operators independent 
of the nucleon velocities: 

Off 1 ' 6 = (1, <n ■ <rj, S tJ ) ® (1, r, ■ Tj ) , (2.4) 

where S^j is the two-nucleon tensor operator. There are only two isospin-conserving spin- 
orbit terms linear in the velocities, with operators: 

Of= 7 > s = L • S ® (1, Ti ■ Tj) , (2.5) 

where L and S are the relative angular momentum and the total spin respectively. The above 
eight terms are unique and able to describe most of the features of the NN interaction. The 
long-range parts of v^Tij) and ^(r^), associated with the <Xj • (Tj Ti ■ Tj and Sij Ti ■ Tj 
operators respectively, are given by the one-pion-exchange potential v n . In addition there 
are phenomenological parameterizations of the short- and intermediate-range parts of the 

Vpijij). 

It is necessary to add several smaller terms to the above eight in order to fit the scattering 
data with a \ 2 ~ 1- These include terms dependent quadratically on the velocity, and static 
and spin-orbit terms breaking the isospin symmetry. In AV18 the quadratic operators are 
chosen as: 



however, in the Paris [19] and Nijmegen [20] models V 2 is used in place of the L 2 



In addition to the isospin-breaking terms in v 7 , strong- interaction isospin-breaking terms 
are necessary to reproduce the data. The Of~ l5 ~ 17 are isotensors with central, <Tj • (Tj, and 
tensor operators, and the long-range part of v p= i5-n(rij) is determined from the difference 
of the masses of charged and neutral pions. The isovector term associated with Of~ ls is 
necessary to explain the difference in the T = 1 , S = pp and nn scattering lengths [Q . The 
interactions associated with the 18 operators listed above contain all the strong and parts 
of the electromagnetic interaction. In addition, there are four more operators that appear 
only in the t> 7 . The number of parameters contained in the AV18 model of is ~ 40, and 
all of them are fairly well determined by the ~ 4300 data in the Nijmegen data base. 
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It is well known that two-nucleon scattering data up to 350 MeV cannot determine the 
potential Vij uniquely. In addition to AV18 , there are four other modern models: Reid-93, 
Nijmegen-I and II |[20|| , and CD-Bonn [2^], all of which fit the Nijmegen data base. The five 
models are different from each other in detail. The Reid-93, Nijmegen-II and AV18 models 
assume that the interaction in each LSJ partial wave can be represented by a local potential 
in that partial wave; in addition the operator structure of the AV18 model given above relates 
the potentials in all partial waves. In states with total spin S — 1 the local potential in LSJ- 
L'SJ coupled waves is expressed as a sum of central, tensor and spin-orbit components. On 
the other hand, the Nijmegen-I and CD-Bonn models include non-local interactions based 
on boson-exchange phenomenology. 

All the models of contain one-pion-exchange potentials, v 71 , as the long-range part, 
and phenomenological shorter-range parts. Fortunately the v n gives the largest contribution 
to nuclear potential energies, and thus the model dependence of the phenomenological parts 
has a limited scope. However, the v n itself is not uniquely predicted by theory. That in the 
CD-Bonn model is derived assuming pseudoscalar pion-nucleon coupling, and is nonlocal, 
while that in the other models is essentially local. The deuteron and 1 So-scattering wave 



functions predicted by the five models are compared in Ref. ||22|| . All the local models 
predict essentially the same wave functions, however, the two nonlocal models give different 
wave functions. The main difference is in the D-state wave function of the deuteron; that 
predicted by CD-Bonn is smaller at r < 2 fm, while that of Nijmegen-I is close to predictions 
of local models at all values of r. The 5-state wave functions of the boson-exchange models, 
CD-Bonn and Nijmegen-I, are larger than those of the local models at r < 1 fm. 

The deuteron elastic-scattering form factors are sensitive to the wave function. The A(q 2 ) 
structure function has been accurately measured, most recently at Jefferson Lab |23|,|24|], 
and results of the recent measurements of the tensor analyzing power, T 2 o(q 2 ), at Jeffer- 
son Lab |25| and NIKHEF |2(J are also available. These indicate that the deuteron wave 
functions calculated from the local potentials are very realistic. They correctly predict the 
observed data with plausible pair currents [fj7| . In addition, it has been shown recently |28 



that nonrelativistic calculations using local v n give deuteron wave functions close to those 
predicted with nonlocal v n obtained with the pseudovector pion-nucleon coupling, favored by 
chiral perturbation theories, and relativistic kinetic energy. The corrections of order p 2 /m 2 , 
to the amplitudes of states with large momentum p, coming from relativistic nonlocalities 
of v n and relativistic kinetic energy cancel in this case. 

The main assumption we make here is that local models provide an accurate represen- 
tation of v^. It is supported by the observed deuteron form factors mentioned above, and 
is valid for the v n as shown in Ref. PS|. The local models also predict essentially the same 
value (7.63 MeV for Reid-93 and 7.62 MeV for AV18 and Nijmegen-II) of the binding en- 
ergy of the triton in non- relativistic calculations with no three-nucleon potential [p9 ]. The 



difference between these values and the observed binding energy of 8.48 MeV is one of the 
indications for the presence of in the Hamiltonian ( |2.1|) . The triton energies obtained 
with the Nijmegen-I and CD-Bonn models, without 8v (Py), are 7.74 and 8.01 MeV respec- 
tively [0. The energies of the alpha particle predicted by AV18 and Nijmegen-II differ by 
only 0.28 MeV while those predicted by Nijmegen-I and CD-Bonn models are respectively 
0.7 and 2.0 MeV more bound than the AV18 value p|. 



Accurate calculations of nuclear matter are not yet practical. Nevertheless the nuclear 
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matter equation of state has been studied for all five modern potentials with the lowest-order 



Brueckner-Hartree-Fock method with continuous single-particle energies ||31|| , again without 
relativistic corrections. The local interactions, Nijmegen-II, Reid-93 and AV18 , give similar 
results, while the most nonlocal CD-Bonn gives the lowest energies. The predicted values 
of equilibrium Eq and po of symmetric nuclear matter are —17.6 MeV at 0.27 fm -3 with 
Nijmegen-II, -18.1 at 0.27 with AV18 , -18.7 at 0.28 with Reid-93, -20.3 at 0.31 with 
Nijmegen-I and —22.9 at 0.37 with CD-Bonn, while the empirical values are —16 MeV at 
0.16 fm -3 . The above Brueckner results for AV18 are quite close to the Eq = —18.2 MeV 
and Pq = 0.3 fm -3 obtained with the variational method using chain summation methods 
The triton and 4 He energies obtained with the nonlocal CD-Bonn interaction are closer to 
experiment than the predictions of local models, but its predicted nuclear matter properties 
are farther away. 



It has been stressed by Friar |32| that the various representations of v n are related 
by unitary transformations. It should be possible to use these transformations to find the 
appropriate current operators that will explain the deuteron form factors with wave functions 
predicted by the nonlocal models. These transformations will also generate three-body forces 
accounting for the difference between energies obtained from local and nonlocal models. 
Thus the deuteron form factors do not exclude nonlocal representations of Vij. However, it 
seems that the simplest realistic models of the nuclear Hamiltonian may be obtained with 
local v^, and fortunately there is much less model dependence in these. In the present 
work we use the AV18 model of v^; however, the other local models will presumably require 
similar V^. 

The two-nucleon interaction depends both on the relative momentum p = (pi — Pj)/2 
and the total momentum P = p; + pj of the interacting nucleons. We can express it as: 

= Vij + 5v(Pij) , (2.7) 

where 5v(P = 0) = 0. The models discussed above give in the P = 0, center of momentum 
frame. In many calculations the is used as an approximation to u y - by neglecting the 
boost correction 5v(Pij). In fact terms dependent on p included in are of the same order 
as those in Sv(Pij) dependent on P [j3~3~| . It is essential to include the 8v(Pij) to obtain 



the true momentum dependence of the For example, the electromagnetic interaction 
between two charges, as well as the analogous vector-meson-exchange interaction between 
two nucleons depends upon px • p 2 = \P 2 — p 2 . The v includes only the p 2 term, while the 
P 2 term is in 8v. The Sv is related to v and its leading term of order P 2 is given by: 

HP) = -^ + -^[P-rP-V^] + -\ [(o-x-^xP-V^] . (2.8) 



The validity of the above equation, obtained by Friar |p4| , in classical and quantum rela- 



tivistic mechanics and in relativistic field theory has been shown in Ref. [33 |. 



The effects of the 5v(Pij) on the energies of 3 H and 4 He [17| and nuclear matter H have 



been studied for the AV18 model using the variational method. This boost correction gives 
a repulsive contribution in both cases. It increases the triton energy by ~ 0.4 MeV away 
from experiment, while the nuclear matter equilibrium Eq and po move to —13.7 MeV at 
0.23 fm -3 , which is closer to the empirical density, but farther from the empirical energy. 



The VMC studies |L7 of 5v (P«) also show that the dominant corrections come from the first 



8 



and second terms of Eq. ( |2.8|) and that only the first six operator terms (the static terms) of 
AV18 give substantial contributions. Accordingly, we ignore the last term of Eq. (|2.8| ) in this 
paper and evaluate the first two for only the static parts of v. Furthermore it was shown 
that the terms arising from the derivatives acting on operators in v were negligible, so we 
do not evaluate them here. 



III. ILLINOIS MODELS OF V ijk 

The Illinois V^u are expressed as: 

V ijk -A 2n U ijk + A 2n U ijk + A 37r U ijk + A R U ijk . [6.L) 

Their four terms represent the V 2n ' PW , V 2n,sw , V 3n ' AR , and V R interactions with strengths 
A™ , A^ , Af R , and A R . In the following subsections we give the spin-isospin and spatial 
operators associated with these interactions and the theoretical estimates of the strengths. 
In the older Urbana models A 1 ™ is denoted by A 2n , A R by U and the V 2n > sw and V 3nAR 
terms are absent. 

A. V 27T ' PW 



The earliest model of V 2n,PW is due to Fujita and Miyazawa |TT|, who assumed that it 



is entirely due to the excitation of the A-resonance as shown in Fig. ^|a. Neglecting the 
nucleon and A kinetic energies we obtain: 

j^PW _ fnNN f^NA m n /o 2) 

2n ~ 81 4tt 4tt (m A -m N ) ' 1 ' ' 

2 ijk W = $Z [{Xij > X jk} {Ti ■ Tj,Tj ■ T k } + ~ [X^, X jk ] [Ti- Tj,Tj ' T k \) , (3.3) 

eye 

= T(m 7r rj i ) Sij + Y(m 1T rij) (Ti ■ <Tj , (3.4) 

Y(x) = — fr(r) , (3.5) 
x 

T(x) = (1 + 1 + 1) Y{x) fr(r) . (3.6) 
\x z x J 

Here £y(r) and £t(?") are short-range cutoff functions. We note that the one-pion-exchange 
two-nucleon interaction used in AV18 is given by: 



2 



l r 

V H = a m - T i ■ T J X n ' ( 3 - 7 ) 



with cutoff functions 



3 4vr 



ty{r) = t T {r) = {\-e-^) , (3.8) 



and c = 2.1 fm -2 . The contact, ^-function part, of the OPEP is not included in the v 77 in 
Urbana- Argonne models since it is difficult to separate it from the other short-range parts. 
These functional forms are used in UIX and all the Illinois models. 



9 



In all the Illinois models except IL3 the cutoff c = 2.1 fm -2 is used and the A™ is varied 
to fit the data, as in UIX. This approximation assumes that the ttAA form factor is similar 
to the ttNN form factor. In fact it is likely that the radius of the A-resonance is larger than 
that of the nucleon, and thus the 7riVA form factor is softer than the irNN. In this case 
use of the T(x) and Y(x) functions from in V 2w,PW would lead to an underestimation 
of A™ ■ In the IL3 model we use a value of A™ typical of the Tucson, Brazil, and Texas 
models and vary the cutoff parameter c in to fit the data. 

Using the observed values of iba and /^va/^ 71 " ~ 0.3, Eq.( |3.2|) predicts that A™ ~ 
-0.04 MeV. With the cutoffs from v ij: the V 2n ' PW of this strength gives a contribution of 
~ —3 MeV to the energy of 3 H. It is much larger than the —0.9 to —0.6 MeV estimated 
by Faddeev calculations |35| , |36[| that include explicit A degrees of freedom and NN ^ AA 
transition potentials |37|]. A part of the difference is probably due to the neglect of kinetic 
energies of the nucleons and A in the energy denominator in Eq.( |3.2|) . Neglecting the 
momenta of the nucleons before the pion emission, the energy denominator in Eq. ( |3.2| ) 
should be wa — + q 2 {l/2m& + l/2mjv), where is the momentum of the first pion in 
Fig. |2|a. The average momenta of pions exchanged in interactions between nucleons in nuclei 
is ~ 500 MeV/c for which Eq. (|3.2| ) underestimates the denominator by ~ 40%. It 

thus appears likely that Eq.( |3.2| ) overestimates the strength of V 2n,PW via the A-resonance 
significantly. 

The other models of Vijk start from the observed pion-nucleon scattering amplitude, and 
use current algebra and PCAC constraints, or chiral symmetry, to extrapolate to the off- 
mass-shell pions responsible for the V 2n . In this way they include the contributions of all 
the irN resonances, as well as that of ttN S-wave scattering to the V 2n . The TM V 2w has 
been cast in the form of Eq. (|3.1|) in Ref . [|lj| . It contains the term V 2n,PW with the operator 
Olj^ FW and the strength A™ = —0.063 MeV. The strength A™ is proportional to the 
parameter b in the ttN scattering amplitude, and the values of b in various models have been 
tabulated by Friar et al. [1141. The Texas model has the largest value of b corresponding to 
A™ = -0.09. These strengths are much larger than —0.04 estimated with the simple Fujita- 
Miyazawa model presumably because of the additional contributions included. However, 
the nucleon and resonance kinetic energies are neglected in the later models, as in the FM, 
therefore their estimate of the strength of V 2n,PW may be too large in magnitude. Another 
concern is that the ttN scattering amplitude used in these models considers only pions of 
momenta less than m n ||14|| . They play a much smaller role in nuclear binding than those 
with momenta ~ 500 MeV/c. 

The factor of 1/4 in the second term of 2 ^ P [Eq. fl3.3|) l, containing the product of 
commutators, is due to the spin and isospin of the A being 3/2. In the TM and later models 
the strength of this term is proportional to the constant d whose values have also been 



tabulated by Friar et al. |14| . The value of d/b is 0.29 in the latest Texas model, however, 



the ratio of the expectation values of the commutator and anticommutator terms of 0^ PW 
is very constant across all the light nuclei studied in this work, and hence this factor cannot 
be determined from the data considered here. We continue to use the Fujita-Miyazawa value 
of 0.25 in this work for simplicity. 
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g y2-K,SW 



The form of the V 2n ' sw , due to ttA S-wave scattering illustrated in Fig. |2|b, in the TM 
model is: 

B(r y , rjfc) {rj ■ Tj, Tj • r fc } {(Sij + <r, ■ crj), (S jk + cr t ■ a k )} , (3.9) 

The BfoijjTjk) contains several terms as given in Ref. [[HJ. We omit the short-range terms 
containing the Z' functions, whose validity has been questioned recently [fTl], and retain only 
the term with pion-exchange-range functions Z[. These functions are given in Eqs.(Al7,A18) 

13 



of Ref. |15|]. The function Z[ is trivially related to the functions Y(x) and T(x) in vf 



[Eqs.( P^| - ^.7D ], and the Z[ contribution to V 2lT ' is expressed as: 



■ SIT / firNN \ i ■-> . .-, , ,1 1 

A 2n = I -^j- J a K , I 3 - 10 ) 

)2n,SW 

ijk 

eye 

Z{x) = | [Y (x) - T(x)} . (3.12) 



2 ijk W = Z(m 7 ,r ij )Z(m 7 ,r jk )(j i ■ r i:j a k ■ r kj Ti ■ r k , (3.11) 



The values of the parameter a' are listed in Ref. ||14|| ; they vary from — 0.51/m n to — 1.87/m^ 
in the recent models. The TM value a' = -1.03/ gives A$j* ~ -0.8 MeV. The V 27r ' sw 
gives rather small contributions to nuclear energies, and it is difficult to extract its strength 
Af^ from nuclear data. In model IL1 we neglect this term, while in all other models it is 
assumed to have the theoretically plausible strength of —1 MeV. 



C. v^^R 

The present model of O^J k AR is based on the three-pion-exchange ring diagrams shown 
in Figs. |2|c,d having only one A at a time in the intermediate states. The V 3n,AR is approxi- 
mated with the sum of V^' R and V^' , which respectively denote the sums of diagrams 
$p and 0d. After neglecting the kinetic energies of the nucleons and the A in the intermediate 
states, diagram |2]c gives: 

Vi,ijk R = Yl / — ; n 2 [vajv-jwW v*U k ) v nn^an{v) + 3 ^ k] , (3.13) 



eye 



lm A - m N 



where vJf N ^ AN (ij), for example, denotes the one-pion-exchange transition potential |37[ 



exciting the nucleon i to the A-resonance state, and j k denotes the term obtained by 
interchanging j and k in the previous term. The above y^ n ' AR can be reduced to a three- 
nucleon operator by eliminating the N ^ A transition spin and isospin operators denoted 
by S, S^, T, and using the generalized Pauli identities: 

T t • T = 2 , (3.14) 
T t xT = -|zr, (3.15) 
T t -AT-B = |A-B-|iT-AxB, (3.16) 
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for the transition isospin operators. The transition spin operators also obey the same iden- 
tities with r replaced with <x. It is useful to reduce the yj* n > AR further by eliminating all 
the terms quadratic in either T\ or a\ (/ = i,j, k) with the Pauli identity: 

<tA(tB = AB + 2(t-AxB, (3.17) 

for <j and r operators. The resulting V^' R contains very many terms which can be 
organized in the following way: 

V%£ R = A&" 0%f k « , (3.18) 

A AR = (\SpL mX |^ 7 - 9 , (3.19) 



Ol*i'jk R ~ 6 [S*S* + A^A 1 ) + 2 ^ (^i^T,ijk + Sl$a,ijk + ^l^aAjk + ^T,ijk^a,ijk) ■ (3.20) 



,ijk \ o r aj ' / < \ a T,ijk ' r a,ijk ' r a,ijk ' r,ijk cr,ijk 

eye 

The letters S and A denote operators that are symmetric and antisymmetric under the 
exchange of j with k. Subscripts r and a label operators containing isospin and spin-space 
parts respectively while superscripts I and D indicate operators that are independent or 
dependent on the cyclic permutation of ijk. The S 1 and A 1 would be more properly written 
with ijk subscripts, but because they are independent of the ordering of ijk, we omit them 
here for brevity. The isospin operators are: 

S T T = 2 + f ( T< • Tj + Tj -T k + T k - Ti) = 4P T=3/2 , (3.21) 

A I T = \ir i -T j ^T k = \ [Ti ■ Tj, Tj ■ Tfe] , (3.22) 

S r,ijk =\ T r T k , (3-23) 

Kiih = > ( 3 - 24 ) 

where we have indicated that is a projector onto isospin 3/2 triples (see the discussion 
at the end of this subsection), and that A T T has the same structure as the commutator part 
of V 27T,PW . The spin-space operators have many terms, and are listed in the Appendix. In 
addition to the spin operators, they contain the functions T(x) and Y(x) in the v*. The 
interaction yj* n < AR has to be symmetric under the exchange of i, j, and k; therefore products 
of S- and A-type operators are not allowed 

The \ 
given by: 



The V^' R , obtained from diagram d in Fig.0, after neglecting the kinetic energies, is 



^2jjk R — X! 7 Y2 i V NA->NN{ik) v AN-*Na(J^) v NN-^NA^j) + J ' ^ ^] > (3.25) 

w (m A - m N ) 2 

After appropriate reductions it can be cast in the form of y^' AR as follows: 

V? AR = A** 0^ R , (3.26) 
f 2 

AAR _ aAR JttNA (r> Q7\ 
JttNN 

n 3Tr,AR _ 8 qlql , 2 Al Al ^STfqlqD ,qIqD i aI aD 1 qD qD \ (QOQ\ 

{ ~ J 2,ijk — 3 D T°a "T z t a 9 \ D a°T,ijk ~t~ J T J a,ijk ~t~ ^r^a^jk 2 T,ijk a,ijk ) ' [o.Zci) 

eye 
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The Vi ' and V^ A are combined using f 2 NA ~ Af 2 ^v^v to obtain: 

^' A * = At* (0%f k R + 40 2 3 ^) = A^O^* , (3.29) 

n 37T,A_R _ 50 ql ql j_ 26 /!/ Al 
^ijk — ~3 D T°a "r ~ /i T /i a 

+ I ^ \^i^?,ijk + ^r^a,ijk + -Ar^ijfc + 13 S^ijkS^jk) • (3.30) 
eye 

The strengths of the terms independent of cyclic permutations are larger than those which 
depend upon them. Therefore we use the simpler V 37T ' AR obtained by neglecting them, i.e. 
with the approximate operator: 

OSf* - fSX + fAiAi . (3.31) 

The value of A AR estimated from the observed values of the constants, and neglecting the 
kinetic energies, is ~ 0.002 MeV. In all the Illinois models the A AR is determined by fitting 
the nuclear energies. 

The V 3n ' AR has an interesting dependence on the total isospin T to t of the three interacting 
nucleons. The can be written as: 

T tot = \ (t< + rj + r k ) , (3.32) 
Si = I T? ot - 1 . (3.33) 

Therefore the first term of V 3n,AR is zero in triplets having T to t = 1/2, i.e., in d + N channels 
as well as in the A = 3,4 bound states. In contrast A ! T is zero in T to t = 3/2 states. It is 
therefore possible to extract the strength of this interaction from the data even though it is 
much weaker than the V 2n . 



D. V 



R 



The pion-exchange three-nucleon interactions are attractive, and lead to significant 
overbinding and large equilibrium density of nuclear matter. Therefore there must be other 
three-nucleon interactions to compensate the attraction from V 2n in nucleon matter at large 
density. In Faddeev calculations of the triton including A excitations [|35,36|, the attraction 
from processes included in V 27T is more than cancelled by "dispersion" terms which describe 
the modification of the contribution of the two-pion-exchange A-box diagrams to due 
to the presence of the third nucleon k. Such repulsive terms also occur in the variational 
theory in which A-excitations are included via transition correlation operators | 40| . The V R 
term in the Urbana models of was designed to approximate these effects. It is retained 
in the Illinois models with the simple spin-isospin independent operator: 



(3.34) 



eye 



35.36 



indicate that the A-effects not included in the 
The value of the strength Ar required to 

0.004 MeV for a 



Results of the Faddeev calculations 
V 2n add ~ 1 MeV to the energy of the triton. 

obtain a contribution of ~ 1 MeV from the V R to the triton energy is 
T(r) with a cutoff c = 2.1 fm -2 . In models IL1-4 the Ar is determined by fitting nuclear 
energies. 



13 



E- V* jk 



Most calculations of light nuclei use the simpler Hamiltonian obtained by approximating 
the Vij in Eq.( |2.1| ) by the Vij. The more accurate Hamiltonian: 

H* = J2 + £(% + MPii)) + E V$ h , (3.35) 

i Zm i i<j i<j<k 

contains the boost correction to the two-nucleon interaction. This correction 8v(Pij) is of 
first order in P?-, and therefore its contribution is calculated as a first-order perturbation in 
wave functions generated by a nonrelativistic Hamiltonian that gives the same final energy 
as H*. In light nuclei the expectation values of 5v and V R are nearly proportional to each 
other [p~6|JT~7[1 . Therefore the energies of light nuclei, calculated with the simpler H can be 
reproduced with the H* using: 

T/* _ APW n 2Tr,PW , A SW n 2Tr,SW , A AR n 3ir,AR , a* n R (Qoa\ 
V ijk — A 2-k U ijk ' U ijk "+* a 3tt U ijk ' ^R^ijk ■ {d.OV) 

which differs from V^k only in the strength of V R . For T{r) with cutoff c = 2.1 fm -2 the 
strength A* R in Vj* k is smaller than Ar by ~ 0.002 MeV. The proportionality of 5v and 
V R contributions appears only in bound light nuclei. In nuclear matter the 5v contribution 
increases more slowly with density than that of V R , and in neutron drops the 5v gives a 
relatively larger contribution. Therefore it is necessary to use the H* in these systems. 



F. The Illinois Model-5 



This model is meant to test the sensitivity of nuclear energies to the spatial shape of V R . 
Here we assume that: 



V$ h = 0.002 R k + A w l[W(i 



(3.37) 



eye 



where W(r) is a modified Woods-Saxon function that has zero- derivative at the origin: 



W{r) 



1 



X _)_ e (r-r w )/aw 



1 + 



r/aw 



1 _|_ e r w la w 



(3.38) 



with rw = 1-0 fm and aw = 0.2 fm. The first term of this V R is meant to take into 
account the contribution of Sv omitted from the simpler H, and the second corresponds to 
the phenomenological three-nucleon repulsion. 



G. Potential Parameters 

The parameters of the Illinois potentials were determined by fitting the observed energies 
of 3 < A < 8 nuclei. The simpler Hamiltonian H was used, and a total of 17 ground or 
excited states with widths less than 200 keV were considered. Table | presents the parameters 
of the new potentials, and for comparison those of UIX Q. The properties of light nuclei 
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calculated from Hamiltonians including the AV18 Vy and the new Vijk are presented below 
in Sec. V. In addition to the various strengths, Table | gives the value, c, of the cutoff 
constant in Eq.(|3.8|) which is used in all the four terms of V^. 

As mentioned earlier, this data set cannot determine all five parameters of the Illinois 
Vijk- At most three parameters were varied for each model, and plausible values are assumed 
for the others. These assumed values are marked with an asterisk in Table |. There is a 
substantial cancellation between the contributions of V 2lT,PW and V R in the nuclear binding 
energies. Therefore one can make correlated changes in the A™ and An, as in models IL2 
and IL4, without significantly spoiling the fit. Presumably nuclear matter calculations with 
these models can help to further constrain the parameters. 

As will be discussed later, the contributions of V 2n,sw and V 2lT ' PW are in a fairly constant 
ratio for the light nuclei considered here, and thus we cannot uniquely determine both A^ 
and A™ . IL1 assumes that A^ = 0, while in all other models A^ = —1 MeV as in 



modern chiral-perturbation-theory potentials |T3|. The A™ in IL2 is less than that in IL1 
by ~ 4% to compensate for the V 2n,sw contribution. All models other than IL3 have the 
same cutoff, c, as in AV18 and UIX and have A™ of the same order of magnitude as the 
Urbana models. This A™ is approximately half that favored by chiral-perturbation theory. 
We constructed the IL3 to see if light nuclei are sensitive to this difference. In this model, 
was fixed at a typical chiral-perturbation theory value and the cutoff parameter, c, 
and strengths A R and A^ p were adjusted to fit the binding energies. The cutoff had to be 
made much softer to compensate for the strong A™ . 

The expectation values of 8v were calculated in a few nuclei for each model, and the 
values of A* R were estimated requiring: 

(Vijk ~ V? jk ) » (5vQ? v )) . (3.39) 

They are listed in Table |. 



IV. QUANTUM MONTE CARLO CALCULATIONS 

The GFMC calculations presented below were made using essentially the same methods 
and variational wave functions as described in Ref. |2| for our calculations of A < 8 nuclei 
with the UIX three-nucleon potential. Here we describe only the few enhancements that 
were necessary for using the Illinois three-nucleon potentials. 

The new terms in the Illinois potentials are static and hence present no formal difficulties 
for the GFMC propagator beyond those already encountered for the UIX V^k- they are 
included by expanding exp(— ^V^Ar) to first order, as in Eq.(4.5) of Ref 0. The structure 
of V 2% ' SW is similar to that of the anticommutator part of V 2n,PW (V 27T,A ) and can also be 
reduced to just two-body operators in spin-isospin space. Thus it can be combined with the 
aV^f in Eq.(4.21) of Ref § with almost no increase in the required computer time. 

The V^ AR involves many operators which are unfortunately not reducible to two-body 
operators in spin-isospin space, and thus, like V 27T,C , the commutator term of V 2n,PW , it adds 
a lot of time to the evaluation of a propagation step. In Ref 0, we showed that V 27T,C could 
be replaced with an increased strength of V 2n ' A in the propagator with no loss of accuracy. 
The V 3lT ' AR cannot be similarly completely replaced, but one can use two to four steps of 
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propagation with just an enhanced V 2n,A and then a correction; see Eq.(4.23) of Ref J2|. In 
general the tests we have made show that systematic errors in the GFMC calculations are 
less than 1% for the total energy; however 8 He appears to be a particularly difficult case 
and the systematic errors for it are probably 2%. 

The GFMC calculations for the models with V^fc and p-shell nuclei were carried out as 
described in Ref. ||. In particular, propagations were made to r = 0.2 MeV -1 with steps 
of At = 0.0005 MeV -1 (400 steps) and expectation values were computed every 20 steps 
with averages of the last 7 values (r > 0.14 MeV -1 ) being used. The s-shell calculations for 
most models and some of the p-shell calculations for models without were propagated 
to only r = 0.1 MeV^ 1 . In most cases 10 unconstrained steps were used before each energy 
evaluation [n u = 10; see Eq.(4.17) of Ref. @], but 20 steps were used for some cases. 

The 4 He energy obtained with the Yakubovsky equations |30| and AV18/UIX Hamil- 
tonian, -28.50(5) MeV, is within ~ 0.5% of the latest GFMC § result, -28.33(2) MeV. 
The GFMC calculations are carried out with an approximation called AV8' [0], containing 
only the first 8 operators given in equations (|2.4|) and fl2.5|) . The small difference between 
the AV18 and AV8' is treated perturbatively. It is therefore likely that the exact results for 
AV18/UIX are a little below the present GFMC results. The differences between GFMC 
and other calculations of binding energies of three and four nucleons will be the subject of 
another paper. 



V. RESULTS - LIGHT NUCLEI 



This section presents various results for 3<v4<8 nuclei and neutron drops using the 
new Illinois models, the AV18/UIX model, and Hamiltonians containing just the AV8' or 
AV18 potentials with no V!y&. All of these results were obtained from GFMC calculations. 
With the exception of the total energy, which is discussed in the next paragraph, the values 
are perturbatively extrapolated from the mixed estimates, as described in Ref. 0. Monte 
Carlo statistical errors are given in parentheses, but no estimate is made for the systematic 
errors associated with the extrapolation of mixed expectation values. The results presented 
in subsections A to F are with the simple H without boost interaction, whose contributions 
are reported in subsection G. 

Most experimental energies and moments are drawn from the standard compilations of 
Ajzenberg-Selove |4l|] and the TUNL Nuclear Data Evaluation Project [42 -Q, while charge 
radii are taken from the NIKHEF compilation ESI. More recent data not included in these 
references are the energy of the 8 He(2 + ) state |46| , the charge radius of 3 He 
A = 8 magnetic and quadrupole moments ESI . 



471, and the 



A. Energies of "Narrow" States 

Table || shows our GFMC energies for 3<A<8 "narrow" states for all the Hamiltonians 
along with the corresponding experimental energies. The states are either particle stable, 
or have experimental widths less than 200 keV, and are used to fit the parameters of the 
Illinois models. The AV8' and AV18 Hamiltonians have no three-nucleon potential; they are 
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presented to show the importance of the and to provide results for comparison with 
those from other many-body methods. 

The AV8' model consists of the Argonne v' s two-nucleon potential, which is an eight- 
operator refit of the AV18 , and the isoscalar Coulomb potential, as defined in Ref. [0. The 
AV8' results should be the most reliable of all our results, because the GFMC propagation is 
made with the same Hamiltonian as is used for the energy expectation values, as discussed 
in Ref. |TJ]. In all other cases an effective H' is used for propagation and a small contribution 
(H — H') is evaluated perturbatively. In the cases that have a three-nucleon potential, the 
Ar is adjusted in H' to make ( H — H' ) ~ 0; however such an ability does not exist for 
the AV18 with no V^fc. The AV18 energies are found by perturbatively evaluating (t> 18 —v' 8 ) 
in the AV8' calculation. 

As described in the previous section, the AV18/IL1 and AV18/IL2 models differ only in 
the values of and A%™, such that ( V 2n > sw + V 2lT ' PW ) in IL2 is nearly the same as 
( Y 2it ,pw ^ j n We made a complete set of calculations for AV18/IL2 and perturbatively 
computed ( V 2n,PW ) for AV18/IL1 in the AV18/IL2 wave functions. This result was used to 
generate 12 of the 21 AV18/IL1 energies in Table [TT]. The procedure was checked by making 
a new propagation using the AV18/IL1 H', and these explicitly calculated values are shown 
for nine states; they are not significantly different from the perturbative estimates, i.e., the 
differences are generally smaller than the Monte Carlo errors. 

The five Illinois Hamiltonians give very similar energies. The predictions for the p- 
shell nuclei are significantly better than those obtained with the AV18/UIX model. This is 
illustrated in Fig. | where results for the AV18/UIX, AV18/IL2, and AV18/IL4 models are 
compared to experiment. The relative stability of the helium and lithium nuclei with the 
Illinois models is clearly evident, as is the just unbound nature of 8 Be. More quantitatively, 
Tables [TIT] and [TV] show various averages of the deviations from experiment for the narrow 
states of Table [0]. Table [TTT| is based on the deviations of the total energies of the 17 states, 



while Table [TV] is based on the deviations of the excitation energies of excited states. Both 
tables show the average deviation (which includes the signs of the deviations), the average 
of the magnitudes of the deviations, and the rms deviations. The average deviations in 
Table [TTT] demonstrate that the Hamiltonians with no systematically underbind these 
nuclei by 5 to 7 MeV; AV18/UIX reduces this to the still large value of 2 MeV underbinding. 
The five Illinois models have no significant systematic under or overbinding. Because the 
errors for the AV8' , AV18 , and AV18/UIX cases are so one-sided, their average absolute 
and rms errors are comparable to their average signed errors. The rms error obtained with 
the AV18/IL1-5 Hamiltonians is ~ 1%. 

Table [TV] shows that the Hamiltonians with no produce a spectrum that is too 
compressed, although it has far smaller deviations than the absolute energies. All of the 
Hamiltonians with Vijk produce excitation spectra with significantly smaller rms deviations. 
However some of them (particularly AV18/IL2) are too expanded. The excitation spectra 
obtained with AV18/IL3 and AV18/IL4 appear to be somewhat better than the others. 

B. Contributions to the Energies 

The contributions of the two-body potentials, Vij, including electromagnetic terms, and 
the sum (K+Vij) are shown for several Hamiltonians in Table [V]. As was discussed in Ref. [|TJ, 
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the perturbatively extrapolated values of (K), (%), and (Vijk) do not add up to the total 
energy, (H). The latter is the most reliably computed quantity. Other studies of GFMC 
calculations suggest that the perturbative extrapolation of the potential energy is more 
reliable than that of the kinetic energy. Therefore the values of (K + v^) are obtained by 
subtracting (Vijk) from (if); estimates of (K) may be obtained by subtracting (vij) from 
(K + Vij ). 

Table [V] shows some of the non-perturbative aspects of these calculations. For the p-shell 
nuclei, the total binding energy steadily increases from AV18 to AV18/UIX to the Illinois 
models. However the values of (K + v^), which involve the same operators in all cases, 
steadily decrease in magnitude. This is because the wave function is being tuned to the 
ever stronger Vijk and hence is becoming less favorable for K + vy. The net increase in the 
binding energy comes from even bigger increases in (Vijk) (s ee Table [VI]) . Although (K + v^) 
is becoming less attractive in this progression, (v^) is becoming more negative due to the 
enhanced tensor-isospin (pion-exchange) correlations induced in the wave function by V^; 
these correlations also increase the kinetic energy. 

Table |V] shows the total (V ijk ) for the various models. The AV18/IL1 and AV18/IL2 
models were constructed to have approximately the same (Vijk)- Although the AV18/IL3 
model has very similar total binding energies as AV18/IL1 and AV18/IL2, there are signifi- 
cant differences in many of its contributions. This indicates that the correlations induced by 
IL3 (which is stronger and has a much softer pion form factor) in the wave functions make 
important non-perturbative changes in (K + Vij). The AV18/IL4 and AV18/IL5 models 
have weaker strengths than the first three and smaller net (Vijk)- As expected, the (Vijk) f° r 
the Illinois models are all larger than for the AV18/UIX model in p-shell nuclei; they are 
also larger for the s-shell nuclei even though all the models give the same binding energies 
for these nuclei. 

The fraction of the total binding energy represented by (Vijk) increases from s- to p-shell 
nuclei and as N — Z increases. For AV18/IL2 it is 17% for 3 H, 30% for 4 He, but then a 
nearly constant 33-37% for 6 ' 7 Li and 8 Be. It then jumps to 49% for 8 Li and 52% for 8 He. 
Expressed as a fraction of the total potential energy, the AV18/IL2 (Vijk) ranges from 2.5% 
for 3 H to 6.1% for 4 He up to 7.5% for 8 He. These fractions are typical of the other Illinois 
models except that AV18/IL3 has somewhat larger ratios. 

The individual contributions to (Vijk) are shown in Tables |V11| and |V111| for AV18/IL2 
and AV18/IL3, respectively. The ratios of the contributions of the anticommutator and com- 
mutator parts of V 2n,PW [see Eq.( |3.3| )] for the nuclei studied here are remarkably constant 
for a given model. For AV18/IL2 the ratio of the anticommutator contribution to the total 
(y2-K,Pw^ var j es only between 0.62 and 0.63 for the nuclei in Table |V11| , while it is 1.0 for 
pure neutron systems (the [r< • Tj, Tj ■ T&] is zero in T — | triples). Very similar ratios are 
obtained for AV18/UIX and the other Illinois models. A similar ratio has also been found in 
VMC calculations of 16 O using an older Urbana model [p0[| . It is because of this very small 
variation that one cannot improve fits to the energies of light nuclei by changing the factor 
of iinEq.p). 

The ratio of the contribution of V 2n,sw to that of V 2lT ' PW is also quite independent of 
nucleus. For AV18/IL2 it ranges from 3.6% for 8 He to 4.2% for 3 H. This ratio depends upon 
the model; for example, it is ~ 2.3% for AV18/IL3 and ~ 5% for AV18/IL4. Thus small 
changes in A™ can accurately compensate large changes in A^ , making it impossible to 
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uniquely determine the value of from the binding energies of light nuclei. We have made 
versions of AV18/IL2 and AV18/IL3 that have the unreasonably large value of = —2.2 
MeV; these also gave good fits to the binding energies of light nuclei. 

The S^S^. and A^A^ terms of V 3n ' AR are shown in the columns labeled V 3n ' ss and 
V 3n ' AA , respectively. The is a projector onto T — | triples and thus the S^S^ term of 
Eq.( |3.3lD vanishes in s-shell nuclei. The A^A^ term results in a repulsive contribution in 
nuclei and hence V 3lT ' AR is repulsive in s-shell nuclei. In p-shell nuclei the S^S^ term of 
y3n,AR j g attractive and larger in magnitude than the A^A^ term. Thus the V 3lT ' AR changes 
sign between s-shell and p-shell nuclei and also becomes more attractive as N — Z increases. 
This allows V 3n ' AR to substantially improve the fit to the energies of light nuclei. Its strength 
is therefore rather well determined by the data in the context of present models. In low- 
density neutron drops the V 2n terms become very small, while the V 3lx ' AR is attractive and 
gives the largest contribution to V^k- 

An interesting property of the A^A^ term is that it only increases by a factor of ~ 3 as 
A increases from 3 to 4, while all the other V^k terms rise by a factor of 5.5. The factor 
of 5.5 can be understood from the simple argument that 4 He has four triples compared to 
only one triple in the trinucleons, and each triple in 4 He gives ~ 40% more contribution 
due to the higher density from increased binding. The V 2n and V R have relatively simple 
spatial dependence given by a product two radial functions; in contrast the A ! a has a more 
complex spatial dependence, given by Eq. (|A6|) , containing radial functions of all the three 
pair distances. In principle, the strength of the A^A^ term in V^k can be adjusted separately 
from the S^S^ term to reproduce the energies of 4 He and the trinucleons with the desired 
precision; however, agreement with less than 1% error can be obtained assuming that the 
strengths of the S^S^ and A^A^ terms of V 3n,AR have the theoretical ratio of 50/26 from 

Table |IX| shows the v n , V 2n , V 3n ' AR , and the remainder contributions to the potential- 
energy expectation values, evaluated for AV18/IL2. Note that additional pion-exchange 
contributions omitted from v n , V 2n and V 3n ' AR are contained in the v R and V R . For example, 
the two-pion-exchange two-nucleon interaction provides most of the intermediate-range part 
of v R . The total two-pion terms of Vijk are typically 11% of the one-pion part of for the 
AV18/UIX, AV18/IL4, and AV18/IL5 models, while they are 16% for AV18/IL2 and 21% 
for AV18/IL3. The ratio of V 3nAR to V 2n chan ges sign between s-shell and p-shell nuclei. 

The ratio of the V 2lT and V R contributions for a given model does not change by much 
in the light nuclei. For the AV18/IL2 model, the ratio is —2.23 for s-shell nuclei and —1.97 
for 8 Li. The AV18/IL4 model has the largest range, from —2.8 for s-shell nuclei to —2.35 
for 8 Li. When only 3 H and 4 He are included in the fit, as in the case of the UIX model, it is 
not possible to determine A™ and Ar separately. For this reason the equilibrium density 
of nuclear matter was used in Ref. [l[ to determine Ar even though exact calculations of 
nuclear matter properties are not yet possible. 

Even after including all A < 8 nuclei in the fit, we cannot determine A™, Ar and A AR 
separately. For example, in AV18/IL2, decreasing Ar and A™ by factors of 0.5 and 0.77 
leaves the energies of 3 H and 4 He unchanged in first order, however, that of 8 Li decreases 
by 0.8 MeV, i.e. by ~ 2%. This change can be compensated by reducing A AR by ~ 20%. 
Thus in first order multiplying A R , A™ and A AR by 0.5, 0.77 and 0.8 leave the energies of 
s-shell nuclei and 8 Li unchanged. 
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The strengths of Ar, A™ and A^J 1 in IL4 are respectively 0.55, 0.76 and 0.81 times 
those in IL2. The difference between these ratios and those in the preceding paragraph 
is due to nonperturbative effects. The overall fit obtained with IL4 is slightly better than 
that with IL2 (see Tables |T| and [TV]), especially for excitation energies, suggesting that the 
parameters of IL4 have more realistic values. Preliminary results obtained for A = 9, 10 
nuclei also indicate that IL4 is better than IL2. 



C. Isobar ic Analog Energies 

Table [X| gives the total isovector (n = 1) and isotensor (n = 2) energy coefficients, a^ T , 
defined in Eq.(5.3) of Ref. ||, for energy differences in isobaric multiplets. The first three 

lines of the table show three different evaluations of the a 1 1 . The first two are expectation 

3,2 

values of the isovector operators in separately computed 3 H and 3 He GFMC wave functions, 
respectively. The last is the difference of separately computed GFMC energies; it has a 
considerably larger statistical error, but otherwise would be the best calculation to compare 
to experiment. The expectation values computed in the two different wave functions are 
statistically different, but, due to its larger error, the energy difference is consistent with 
both. 

The o^t m multiplets having A > 3 (Table |X|) have been computed as expectation 
values of the isovector and isotensor parts of the Hamiltonians in the wave function for 
the Z = ^A — T nucleus. Table [Xl| gives a breakdown of these coefficients into various 



contributions for the AV18/IL2 model. These include: K CSB , the kinetic energy contribution 
due to the neutron- proton mass difference; Vci(pp), the proton-proton Coulomb term; v 1,R , 
the remaining electromagnetic contributions, such as magnetic moment interactions, which 
are part of Argonne v\g] and the strong-interaction terms v CSB and v CD . The isotensor v CD 
comes from components 15—17 of Argonne V\s and contributes only to the a^ 2 \ while the 
isovector v CSB is term 18 and contributes only to a^. 

The isovector terms are dominated by vci(pp), the expectation value of which is strongly 
correlated with the rms radii. The Illinois models generally give better total binding energies 
and radii (see below) and thus better values for these coefficients. However, the remaining 
kinetic and potential terms contribute 5-10% of the total isovector energy coefficient, thus 
playing an important part in the overall good agreement with experiment. The isotensor 
terms are also dominated by vci(pp), but in this case the increased binding of the Illinois 
models has not improved the agreement with experiment. The vci{pp) alone underestimates 

(2) 

the a\ T ; the strong interaction contributions have the correct sign, but seem to be too large 
in magnitude. 



D. Energies of "Wide" States and Spin-orbit Splittings 

So far we have presented results for states of nuclei that are either particle stable or 
have narrow experimental widths. Our GFMC calculations which treat all states as bound 
systems should be reliable for such states, and the comparison of the resulting energies 
to experimental values should be unambiguous. Table [XII| shows GFMC calculations of 
additional states that are experimentally broad or experimentally unknown. An example 
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of the GFMC propagation for broad and narrow states was shown in 0; the energy of the 
broad state falls slowly but steadily with imaginary time. This introduces a certain, but 
usually small, ambiguity in determining the resonance energy from the calculation; we use 
the average of the the energies for 0.14< r <0.20 MeV -1 . In addition, the experimental 
assignment of some resonance energies may be difficult. Nonetheless, the rms errors in the 
Illinois predictions of the experimentally known states in Table [XII| are only ~700 keV. 

It has long been known that Hamiltonians containing only realistic two-nucleon potentials 
often cannot correctly predict the observed spin-orbit splitting of nuclear levels; in fact one 
of the original motivations for the Fujita-Miyazawa three-nucleon potential was the study of 
spin-orbit splittings fllf . In [[51]] we showed that one of the Urbana family of makes a 



substantial contribution to the spin-orbit splitting in 15 N. Table |XII1] shows calculated and 
experimental splittings for a number of states that are spin-orbit partners in conventional 
shell-model calculations. The dominant L and S in the shell-model calculations (and in the 
one-body parts of our variational wave functions) are shown. 

The spin-orbit splitting computed with just the two-nucleon interactions, AV8' or 
AV18 are generally too small, sometimes by factors of two to three. In some cases AV18/UIX 
makes a significant increase, but in general it also predicts too small splittings. The predic- 
tions with the Illinois models are in much better agreement with the experimental values. 
Due to significant statistical errors in the calculated spin-orbit splittings and the fact that 
some of the spin-orbit partners are wide states, they cannot yet be used to differentiate 
between the various Illinois models. 



E. Point Nucleon Radii and Electromagnetic Moments 

Table [XI V| gives the point proton and neutron radii for some of these models. The 
"experimental" point proton radii were obtained by subtracting a proton mean-square radius 
of 0.743 fm 2 and times a neutron mean-square radius of -0.116 fm 2 from the squares of 
the measured charge radii. As mentioned earlier the GFMC propagations are carried out 
for an H', and results for the desired H are obtained by treating H — H' as a first order 
perturbation. Therefore the point radii and electromagnetic moments are computed for H' 
instead of H . For the models with three-body potentials, the H' has a modified An such 
that ( H — H' ) ~ 0, and thus the radii and moments for the H' should be close to the 
desired ones for H. However the calculations for AV18 with no three-body potential use 
just AV8' for H'. In this case ( H — H' ) is significantly different from zero, and we can 
only quote radii and moments for AV8' . For the few light p-shell nuclei that have measured 
radii, those obtained with the Illinois models, which produce good binding energies, are in 
better agreement with the data than those obtained with either AV8' or AV18/UIX. 



Table KV| shows the experimental isoscalar and isovector magnetic moments for the cases 
that have been measured along with values calculated using only one-body current operators. 
The values in the table are defined in the same way as the coefficients a^ T in Table [X] (see 
Eq.(5.3) of Ref. ||) and thus the isovector (/x^) values for T = | cases are —2 times those 
often quoted. The fj,^ and /i^ are obtained from expectation values of the isoscalar and 
isovector magnetic-moment operators in the wave function for the nucleus having smallest 
Z = |v4 — T. In this approximation the isotensor = 0, since the one-body magnetic- 
moment operator does not have an isotensor term. However, one may obtain a small /i*- 2 ** 
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when the magnetic moments are separately calculated for each state in the multiplet due to 
violation of isospin symmetry in the wave functions. 

For the A = 8, T = 1 nuclei, the experimental isoscalar and isovector moments are 
obtained from the sum and difference of the values for B and Li, since the magnetic moment 
of the T = 1 J T = 2 + state in Be is not measured. In fact the sum gives 2//°) + fi^ 2 \ but in 
the present approximation /i^ = 0. 

The computed magnetic moments show little dependence on the three-nucleon inter- 
action. Because pair currents are not included in the calculated values in Table [XV] , one 
cannot expect good agreement with the experimental values, especially for the isovector 
values. The pair-current corrections computed for the A = 3 system using the AV18/UIX 
model are 0.034 and —0.778 for the isoscalar and isovector moments, respectively f52[. This 
isoscalar correction is twice what is needed to achieve agreement with experiment while the 
isovector value results in perfect agreement. All the computed magnetic moments differ from 
the experimental values by amounts comparable to the corrections computed for A = 3. 

Table [XVI] shows computed (again using just impulse approximation) and experimental 
quadrupole moments. The Illinois models predict quadrupole moments that are generally 
smaller than those obtained using AV18/UIX. This is a consequence of the increased binding 
energy, and resulting smaller rms radii. The situation for 6 Li is very difficult due to the large 
cancellation between orbital and intrinsic components (in an alpha-deuteron picture) of the 
quadrupole moment. 



F. Neutron Drops 

Neutron drops are systems of interacting neutrons confined in an artificial external well. 
We have previously reported results for systems of seven and eight neutrons as a basis for 
comparing Skyrme models of neutron-rich systems with microscopic calculations based on 



realistic interactions |L8[ . The determination of the isospin dependence of the Skyrme model 
spin-orbit parameters is of particular interest. The external one-body well that we use is a 
Woods-Saxon: 

^>^ i + e* Pl 4U)/ ao r (51) 

the parameters are Vo = —20 MeV, r = 3.0 fm, and = 0.65 fm. Neither the external 
well nor the total internal potential (%■ + V^) are individually attractive enough to produce 
bound states of seven or eight neutrons; however the combination does produce binding. 

Many of the tables show results for the neutron drops. The T — | nature of the S^S^. term 
of v 3n,AR results in large contributions in the neutron drops. As a result the seven-neutron 
drops computed with some of the Illinois potentials have double the spin-orbit splitting 
predicted by AV18/UIX. This strong dependence on the Hamiltonian indicates that the 
conclusion of [|L8|| , that conventional Skyrme models over-predict the spin-orbit splitting in 



neutron-rich systems may not be valid. The IL3 model gives larger spin-orbit splitting in 
the A = 7 neutron drop than the others. 
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G. 5v Contributions 



Table [XVI H shows the expectation values of Sv in various nuclei for some of the Illinois 



Hamiltonians along with the net change in the binding energy due to the boost correction: 

(H*-H) = (6v + (V* jk -V ijk )) . (5.2) 

In the light nuclei the net change is, at most, comparable to 1% of the binding energies, and 
therefore the Hamiltonian H* [Eq.( |3.35" )1 gives essentially the same energies as the simpler 



H [Eq.( [2.1|) l without 5v correction. However, the net change is not necessarily small in 
other nuclear systems like nuclear or neutron matter or neutron drops. In the A = 8 drop 
(Table |XVII| ) the (V^k — V*j k ) is only half as large as (6v). In these systems we must use 



the Hamiltonian H* 



VI. CONCLUSIONS AND DISCUSSION 

We find that nonrelativistic Hamiltonians containing two- and three-nucleon interactions 
can reproduce the energies of all the bound and narrow states of up to eight nucleons with 
an rms error < 1% via GFMC calculations, which have an estimated error of < 2%. The 
three-nucleon interactions give a significant fraction of nuclear binding energy due to a large 
cancellation between the kinetic and two-nucleon interaction energies. This cancellation is 
seen even in the deuteron whose kinetic and interaction energies obtained from AV18 are 
respectively +19.9 and —22.1 MeV. Since the AV18 model is very successful in explaining 
the observed deuteron form factors, and all realistic models of also have this feature, it 
seems to be inherent in nuclei. 

The Bochum group |53[ obtained a one parameter family of Vijk models by choosing the 



cutoff mass in the TM model of V 2n , including P- and S-wave terms, to reproduce the 3 H 
or 3 He energy with several modern «y. More recently they have obtained results for 



3 H, 3 He and 4 He with AV18 and a revised version of TM V 27T as well as the UIX model. 
They observe that when the energy of 3 He is reproduced that of 4 He is very close to the 
observed, but 3 H is underbound by ~ 0.5%. On the other hand if 3 H energy is reproduced 
both 3 He and 4 He are overbound by less than 1%. In order to improve the accuracy of 
nuclear Hamiltonians beyond 99%, a more quantitative description of the charge symmetry 
breaking interactions is necessary. 

The above results limit the contribution of four-nucleon interactions in 4 He to less than 
1% of its binding energy. The (Vijk) contributes approximately a quarter of the 4 He binding, 
although the difference of the total energies computed with and without a V^k is only a 
seventh of the binding energy. The ratio of the expectation values of V^k and Vij is ~ 7% in 
the A = 8 nuclei, and (V^k) gives up to one half of their binding energy. If the four-nucleon 
interactions were to contribute ~ 6% of V^k, they could influence the A = 8 binding energies 
by ~ 3%. The present calculations do not indicate a need to include four-nucleon interactions 
to fit the observed energies at the 1% level. Thus, either the four-nucleon contributions are 
smaller than 1% of the binding energies, or parts of the present models of V^k are mocking 
up their effects. 

The energies of light nuclei can be used to determine at most three parameters of V^k- 
We can choose them as either the strengths A™ , A^ and Ar, or use a theoretical value 
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of A™ an d fit the short-range cutoff. It is possible to make correlated changes in the three 
strengths, as in models IL2 and IL4, which have relatively small effect on the energies of 
light nuclei. 

All the realistic models of V^, including the older Urbana models, have a cancella- 
tion between the attractive V 2lT and repulsive V R . For example, their contributions in the 
AV18/IL2 model, to the energy of 3 He ( 8 Be) are respectively -2.8 (-37) and 1.3 (18) MeV. 
The contribution of V R grows faster than that of V 2n as either A or the density of matter 
increases, and lowers the saturation density of symmetric nuclear matter. It is difficult to 
determine the Ar quantitatively from the energies of light nuclei. 

Such a cancellation was noticed in Faddeev calculations including A components in 



the triton wave function f35U3(| . The estimates of Picklesimer, Rice and Brandenburg |36 
indicated that the A part of V 2 ' 7T,PW changes the 3 H energy by only ~ —0.75 MeV, while 
the processes, which represent the suppression of the attractive two-nucleon v 2n by the third 
nucleon via the dispersion effect, give ~ 1.1 MeV. The V R contribution in present models 
is comparable to their A dispersive effect, while the V 2lT ' PW is much more attractive. 



Studies of A = 3, 4 nuclei with relativistic Hamiltonians |I7j indicate that the boost 
interaction 5v gives the largest relativistic correction, of ~ 0.4 MeV, to the triton energy. 
It is included in the present relativistic {H*) models. The other corrections included in 
the relativistic Hamiltonians, but excluded here, are only ~ 0.1 ± 0.05 MeV, i.e. of order 
1%. However, three-nucleon interactions via Z-diagrams fl54 |, if any, have to be added to 



the V*j k in the relativistic Hamiltonians. Forest |55| has estimated their contribution to 



the triton, using the scalar and vector parts of AV18 Vij, obtained with Riska's method, 
to be ~ 0.3 MeV. In the present Illinois models these are also buried in the V R *. In the 



initial Illinois models discussed at the International Nuclear Physics Conference in Paris |56 
we attempted to include possible Z-diagram contributions in Vijk, however, the observed 
energies can be reproduced without them. If these exist, then, as in the case of V 27T,SW , we 
can assume theoretical values for their strengths and fit the energies, presumably with ~ 1% 
accuracy, by readjusting the strengths A™, A^ R and A* R . Nuclear binding energies seem 
to require only three components in the V*j k , an attractive part to provide more binding to 
light nuclei, a repulsive part to make nuclear matter E(p) saturate at empirical po, and an 
isospin-dependent term to provide extra binding to the neutron-rich helium isotopes. 

Preliminary versions of the IL1-IL3 models were also discussed at the Few-Body Physics 
conference in Taipei |57|. The parameters of the Illinois models reported there were incorrect 



due to a programming error. The correct values are as reported here. Since then the IL3 
model has been revised. In the present IL3 model, softer cutoffs are used also in the V R 
terms. Results for nuclear and neutron matter with the Illinois models will be published 
separately. 

As mentioned in the introduction, additional data necessary to further the study of V^k 
may be obtained from the scattering of nucleons by deuterium. It is known that most of 
the low energy Nd elastic scattering observables are well reproduced by realistic models of 

WJ\. The expectation value of is only ~ 2.5% of that of + Vjk + v ki in 3 H, thus it 
is not expected to have a large effect on this scattering. However, all realistic models of 
underestimate the observed nucleon analyzing power A y in low energy Nd scattering; the 
spin-orbit splitting induced by the V 2n ' PW of the various Vijk reduces the error somewhat 
but the additional spin-orbit splitting induced by the present Illinois models is probably 
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inadequate to completely correct the underestimate |58[ . 

More recently it has also been suggested that minima and polarization observables of Nd 
elastic scattering at intermediate energies are sensitive to Vyk, and may be used to refine 
models of V^k [ 59| — |61f1 . Thus further improvements in realistic models of V^k may be possible 
by a simultaneous fit to Nd scattering observables and nuclear binding energies. Since Nd 
scattering is sensitive only to the Vijk in the total isospin T = 1/2 state, it does not provide 
information on the Vijk in the T = 3/2 state. This work shows the need for an attractive 
interaction in T = 3/2 states to reproduce the energies of p-shell nuclei. It may be possible 
to access this part of V^k in n — t and p— 3 He scattering. 
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APPENDIX: THE SPIN-SPACE OPERATORS IN V 37r ' AR 

It is convenient to define the functions: 

3T(m n rij) 



t 



ij 



yij = Y{m % rij) - T(m lv r ij ) , 
with which the spin-space operator Xy in vfj [Eq.( |3.7| )1 can be expressed as: 

ij tij&i ' ^-ij & j ' ^ij ~t~ VijC i ' & j • 

Evaluating the A-ring diagrams using closure approximation then gives: 

S a = ^yijVjkVki + | ^^( r ijtijyjkyki + C ' jtijtjkDki) ~ |C \C ' jC 'ktijtjktki 



eye 



+ 



E 



(T : ■ (J a 



eye 



\ /jGj ' O" kCjtijtjkVki 



eye 



(Al) 
(A2) 

(A3) 



yijVjkVki + | J2 Aj^iiVikVu 

eye 

3 ^ ^ (""j ' Tjj (Tj • T{j tijUkiU jk &i ' I"fci &j ' Tfcj tkiVjkVij &i ' ^ jk & j ' ^jk tjkVijVki) 
eye 

3 ^ t Cj^CTi ■ TTjfc (Tj • YkitkitjkVij 
eye 

+ \ E] <T ' ' a ' a {tijtjkVki + tijUjktki + Cktijtjktki) , (A4) 



eye 
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S a,ijk = \ a 3 ' a k [ZyijVjkyki + CfUji/jktki + rfjtijykiyjk + r 2 jk t jk yijy ki + r^^-fc^ 

~\~ ~jf-^i & 3 ' ^ki &k ' ^ij t'ijy jktki 3^"j ' ^ij ^k ' ~^ij tijykiyjk 
3^3 ' ""^ ' r ^ ^kiVjkVij s^j ' 1 ' jk & k ' ^ jk tjkVijVki 

3°"i ' a ""^ ' a {tijtjkyki + yijtjktki + C{tijt jktki) j (A5) 

= ^[cr i • cr j x cr k y^y^Vki + cr i acTjacr k a tijt jk t ki ] 

+ (^i X ^ j ' r ij ^k ■ ?ij UjVjkVki + (Ti ■ a (Tj ■ <T k Citijy jk t ki ) 

eye 

+ 3 ' T 3 k <Tk ' r *i ^3 ' a tijtjkVki 

eye 

+ y ^ " a (Citijyj k t k i — Cjtijtj k y k i — C k yijtj k t k i — CjC k tijtj k t k j) , (A6) 

eye 

and 

^<r,ijk = ~3 (T i ' a (CitijV jktki ~ Cjtijtj k yki ~ C k yijt jktki ~ C ' jC 'ktijt jktki) ■ (A7) 

The Ci = Tij ■ r ik and a = x Vj k . A cyclic sum over the indices ijk is denoted by eye. 
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TABLES 



TABLE I. Three-body potential parameters used in this paper. Parameters that were not 



varied in 


fitting 


the data are 


marked with 


an *. 








Model 


c 


aPW 


aSW 

A 2n 


aAR 
A 3n 


A 

A R 


A 

A w 


A * 

A R 




fm~ 2 


MeV 


MeV 


MeV 


MeV 


MeV 


MeV 


UIX 


2.1* 


-0.0293 






0.00480 


0.* 


0.00291 


ILl 


2.1* 


-0.0385 


0.0* 


0.0026* 


0.00705* 


0.* 


0.00491 


IL2 


2.1* 


-0.037 


-1.0* 


0.0026 


0.00705 


0.* 


0.00493 


IL3 


1.5 


-0.07* 


-1.0* 


0.0065 


0.032 


0.* 


0.02562 


IL4 


2.1* 


-0.028* 


-1.0* 


0.0021 


0.0039 


0.* 


0.00196 


IL5 


2.1* 


-0.03 


-1.0* 


0.0021* 


0.002* 


210 


0.0 



TABLE II. Experimental and GFMC energies (in MeV) of particle-stable or narrow-width 
nuclear states and of neutron drops. Monte Carlo statistical errors in the last digits are shown in 
parentheses. The final column gives experimental widths in keV. 





AV8' 


AV18 


UIX 


ILl 


IL2 


IL3 


IL4 


IL5 


Expt. 


r 


3 H(i + ) 


-7.76(1) 


-7.61(1) 


-8.46(1) 


-8.43(1) 


-8.43(1) 


-8.41(1) 


-8.44(1) 


-8.41(1) 


-8.48 




3 He(i + ) 


-7.02(1) 


-6.87(1) 


-7.71(1) 


-7.68(1) 


-7.67(1) 


-7.66(1) 


-7.69(1) 


-7.66(1) 


-7.72 




4 He(0+) 


-25.14(2) 


-24.07(4) 


-28.33(2) 


-28.38(2) 


-28.37(3) 


-28.24(3) 


-28.35(2) 


-28.23(2) 


-28.30 




6 He(0+) 


-25.20(6) 


-23.9(1) 


-28.1(1) 


-29.4(1) 


-29.4(1) 


-29.3(2) 


-29.3(1) 


-29.5(1) 


-29.27 




6 He(2+) 


-23.18(6) 


-21.8(1) 


-26.3(1) 


-27.2(1) 


-27.1(1) 


-27.8(1) 


-27.4(1) 


-27.3(1) 


-27.47 


113 


6 Li(l+) 


-28.19(5) 


-26.9(1) 


-31.1(1) 


-31.9(1) 


-32.3(1) 


-32.2(1) 


-32.0(1) 


-32.1(1) 


-31.99 




6 Li(3+) 


-24.98(5) 


-23.5(1) 


-28.1(1) 


-30.1(2) 


-30.1(2) 


-30.0(2) 


-29.8(2) 


-29.8(2) 


-29.80 


24 


7 He(f-) 


-22.82(10) 


-21.2(2) 


-25.8(2) 


-29.3(3) 


-29.2(3) 


-29.3(3) 


-29.3(3) 


-29.2(2) 


-28.82 


160 


7 Li(|") 


-33.56(6) 


-31.6(1) 


-37.8(1) 


-39.4(2) 


-39.6(2) 


-39.3(2) 


-39.5(2) 


-39.3(2) 


-39.24 






-33.17(7) 


-31.1(2) 


-37.5(2) 


-39.2(2) 


-39.1(2) 


-38.7(2) 


-39.0(2) 


-39.0(2) 


-38.77 




'Li(f-) 


-28.41(6) 


-26.4(1) 


-32.1(1) 


-34.5(3) 


-34.4(3) 


-34.0(2) 


-34.5(2) 


-34.2(3) 


-34.61 


93 


8 He(0+) 


-23.8(1) 


-21.6(2) 


-27.2(2) 


-30.5(3) 


-31.3(3) 


-32.0(4) 


-31.9(4) 


-31.0(2) 


-31.41 




8 Li(2+) 


-34.2(1) 


-31.8(3) 


-38.0(2) 


-41.8(3) 


-42.2(2) 


-41.2(3) 


-42.0(3) 


-42.5(3) 


-41.28 




8 Li(l+) 


-33.9(1) 


-31.6(2) 


-37.4(2) 


-40.5(3) 


-40.5(3) 


-40.2(3) 


-40.9(3) 


-40.9(3) 


-40.30 




8 Li(3+) 


-31.4(1) 


-28.9(2) 


-35.3(2) 


-39.3(3) 


-39.1(3) 


-39.5(4) 


-39.3(3) 


-39.2(3) 


-39.02 


33 


8 Li(4+) 


-28.1(1) 


-25.5(2) 


-31.7(2) 


-34.9(3) 


-35.0(3) 


-34.7(3) 


-35.2(3) 


-34.9(3) 


-34.75 


35 


8 Be(0+) 


-47.9(1) 


-45.6(3) 


-54.4(2) 


-57.2(4) 


-56.6(4) 


-55.6(4) 


-56.5(3) 


-55.7(3) 


-56.50 




8 Bc(l+) 


-32.8(2) 


-30.9(3) 


-36.3(3) 


-37.8(2) 


-37.6(2) 


-37.3(3) 


-38.8(3) 


-38.9(4) 


-38.35 


138 


7 n(D 


-33.78(4) 


-33.47(5) 


-33.2(1) 


-36.0(2) 


-35.8(2) 


-36.6(3) 


-35.2(3) 


-35.3(3) 






7 n(|") 


-32.25(4) 


-31.82(5) 


-31.7(1) 


-33.2(2) 


-33.0(2) 


-33.0(3) 


-32.9(3) 


-33.1(2) 






8 n(0+) 


-39.73(6) 


-39.21(8) 


-37.8(1) 


-41.3(3) 


-41.1(3) 


-40.7(2) 


-40.7(3) 


-40.7(2) 
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TABLE III. Average deviations (in MeV) from experimental energies. For each Hamiltonian, 
the average signed deviation, average magnitude of deviation and rms deviation are shown for the 
17 "narrow" states given in Table |n] (Only 3 He energies are used for A = 3). 



Model 


Average 


Average 


rms 




Deviation 


Deviation 


Deviation 


AV8' 


5.52(2) 


5.52 


5.83 


AV18 


7.32(5) 


7.32 


7.72 


AV18/UIX 


2.02(4) 


2.02 


2.34 


AV18/IL1 


-0.09(6) 


0.31 


0.38 


AV18/IL2 


-0.10(6) 


0.28 


0.36 


AV18/IL3 


0.04(7) 


0.31 


0.44 


AV18/IL4 


-0.21(6) 


0.24 


0.33 


AV18/IL5 


-0.12(6) 


0.34 


0.46 



TABLE IV. Average deviations (in MeV) from experimental excitation energies for the 8 "nar- 
row" excited states. As in Table ITJ, but for excitation energies rather than total energies. 



Model 


Average 


Average 


rms 




Deviation 


Deviation 


Deviation 


AV8' 


-0.23(5) 


0.83 


1.20 


AV18 


-0.22(10) 


0.90 


1.36 


AV18/UIX 


0.17(8) 


0.41 


0.53 


AV18/IL1 


0.29(13) 


0.44 


0.53 


AV18/IL2 


0.53(12) 


0.53 


0.61 


AV18/IL3 


0.03(14) 


0.24 


0.34 


AV18/IL4 


0.09(12) 


0.20 


0.25 


AV18/IL5 


0.27(13) 


0.66 


0.79 
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TABLE V. Ground-state expectation values of the two-body potential and (K + Vij) (in MeV) 
for the AV8' , AV18 , AV18/UIX, and AV18/IL2 Hamiltonians. 





AV8' 


AV18 




AV18/UIX 


AV18/IL2 




Vij 


Vij 


K + 


Vij 


K + 


Vij 


K + v^ 


3 H 


-54.9(2) 


-54.8(2) 


-7.605(5) 


-58.7(2) 


-7.27(1) 


-58.6(2) 


-6.97(1) 


3 He 


-53.4(2) 


-53.3(2) 


-6.882(5) 


-57.1(2) 


-6.54(1) 


-56.7(2) 


-6.27(1) 


4 He 


-126.(1) 


-124.9(7) 


-24.07(4) 


-135.9(5) 


-21.98(6) 


-136.4(5) 


-19.99(8) 


6 He 


-155.(1) 


-153.(1) 


-23.89(8) 


-164.(1) 


-21.1(2) 


-171.(2) 


-17.9(3) 


6 Li 


-174.(1) 


-173.(1) 


-26.9(1) 


-182.(1) 


-23.8(2) 


-187.(2) 


-21.2(3) 


7 Li 


-221.(2) 


-219.(2) 


-31.6(1) 


-225.(2) 


-28.7(2) 


-232.(3) 


-25.1(5) 


8 He 


-193.(3) 


-191.(3) 


-21.6(2) 


-194.(1) 


-19.1(2) 


-218.(3) 


-15.(1) 


8 Li 


-249.(4) 


-247.(3) 


-31.8(3) 


-255.(2) 


-27.8(3) 


-278.(2) 


-21.6(4) 


8 Be 


-287.(3) 


-284.(3) 


-45.6(3) 


-297.(2) 


-39.6(4) 


-303.(3) 


-35.5(8) 


7 n 


-54.8(7) 


-54.5(7) 


-33.47(5) 


-54.2(8) 


-33.86(9) 


-59.(1) 


-32.2(4) 


8 n 


-69.8(6) 


-69.3(6) 


-39.21(8) 


-65.9(8) 


-38.8(1) 


-73.(1) 


-38.2(5) 



TABLE VI. Total three-nucleon potential energies (in MeV) for the AV18/UIX and Illinois 
Hamiltonians. 





AV18/UIX 


AV18/IL1 


AV18/IL2 


AV18/IL3 


AV18/IL4 


AV18/IL5 


3 H 


-1.19(1) 


-1.46(1) 


-1.46(1) 


-1.65(1) 


-1.25(1) 


-1.24(1) 


3 He 


-1.17(1) 


-1.44(1) 


-1.41(1) 


-1.64(1) 


-1.22(1) 


-1.24(1) 


4 He 


-6.35(5) 


-8.4(1) 


-8.38(7) 


-10.02(7) 


-7.18(6) 


-7.24(5) 


6 He 


-7.0(1) 


-11.3(2) 


-11.5(3) 


-14.1(3) 


-9.8(2) 


-9.9(2) 


6 Li 


-7.3(2) 


-11.2(3) 


-11.1(3) 


-13.6(3) 


-9.6(2) 


-10.2(2) 


7 Li 


-9.1(2) 


-13.8(4) 


-14.5(4) 


-16.9(4) 


-12.8(4) 


-13.4(4) 


8 He 


-8.0(2) 


-15.5(5) 


-16.3(5) 


-19.7(6) 


-16.(1) 


-14.7(3) 


8 Li 


-10.2(2) 


-19.3(5) 


-20.6(4) 


-24.9(6) 


-18.2(4) 


-17.3(5) 


8 Be 


-14.9(3) 


-20.8(7) 


-21.(1) 


-25.1(8) 


-19.0(4) 


-19.4(5) 


7 n 


0.69(4) 


-3.8(3) 


-3.6(3) 


-4.8(4) 


-2.5(3) 


-2.5(4) 


8 n 


1.01(6) 


-3.1(4) 


-3.0(4) 


-3.7(3) 


-2.8(3) 


-1.9(2) 
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TABLE VII. Contributions of various three-body potential terms (in MeV) evaluated for the 
AV18/IL2 Hamiltonian. 





y2n,A 


y2n,C 


y2n,SW 


y3n,SS 


y3n,AA 


yR 


3 H 


-1.78(1) 


-1.082(9) 


-0.119(1) 





0.182(3) 


1.34(1) 


3 He 


-1.72(1) 


-1.045(8) 


-0.115(1) 





0.176(3) 


1.29(1) 


4 He 


-9.76(8) 


-5.85(5) 


-0.652(5) 





0.63(1) 


7.26(7) 


6 He 


-12.2(3) 


-7.3(1) 


-0.74(2) 


-1.33(6) 


0.42(4) 


9.6(3) 


6 Li 


-11.9(2) 


-7.2(1) 


-0.72(2) 


-0.81(4) 


0.37(4) 


9.1(2) 


7 Li 


-15.4(4) 


-9.3(2) 


-0.91(3) 


-1.61(9) 


0.5(1) 


12.3(4) 


8 He 


-15.6(4) 


-9.1(2) 


-0.88(3) 


-4.5(2) 


0.44(7) 


13.3(4) 


8 Li 


-20.7(3) 


-12.2(2) 


-1.26(2) 


-4.4(1) 


0.55(5) 


17.4(3) 


8 Be 


-23.1(6) 


-14.0(3) 


-1.38(4) 


-1.6(1) 


0.7(1) 


18.3(6) 


7 n 


-0.15(5) 




0.08(1) 


-5.4(3) 





1.9(1) 


8 n 


0.13(9) 




0.18(1) 


-5.9(4) 





2.6(2) 



TABLE VIII. Contributions of various three-body potential terms (in MeV) evaluated for the 
AV18/IL3 Hamiltonian. 





y2w,A 


y2n,C 


y2n,SW 


y3Tr,SS 


y3n,AA 


yR 


3 H 


-2.37(2) 


-1.47(1) 


-0.092(1) 





0.217(3) 


2.07(2) 


3 He 


-2.35(2) 


-1.46(1) 


-0.091(1) 





0.212(3) 


2.04(2) 


4 He 


-12.99(9) 


-7.95(5) 


-0.525(4) 





0.54(1) 


10.9(1) 


6 He 


-16.4(3) 


-9.9(2) 


-0.61(1) 


-1.9(1) 


0.26(4) 


14.5(3) 


6 Li 


-16.2(3) 


-9.9(2) 


-0.59(1) 


-1.31(6) 


0.25(5) 


14.2(4) 


7 Li 


-20.9(5) 


-12.8(3) 


-0.73(2) 


-2.32(9) 


0.38(7) 


19.5(6) 


8 He 


-21.6(5) 


-12.6(3) 


-0.78(3) 


-7.0(2) 


0.16(9) 


22.2(7) 


8 Li 


-27.8(6) 


-16.4(3) 


-0.96(2) 


-6.2(2) 


0.13(9) 


26.3(7) 


8 Be 


-30.8(7) 


-18.9(4) 


-1.12(3) 


-2.4(2) 


0.7(1) 


27.4(9) 


7 n 


-0.13(8) 




0.07(1) 


-8.5(5) 





3.8(2) 


8 n 


0.52(8) 




0.17(1) 


-8.9(4) 





4.5(1) 
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TABLE IX. Contributions of two-nucleon and three-nucleon pion and remainder potentials (in 
MeV) for the AV18/IL2 Hamiltonian. 





V 7T 


V R 


V 2 * 


y3ir 


yR 


3 H 


-45.0(2) 


-13.5(2) 


-2.98(2) 


0.182(3) 


1.34(1) 


3 He 


-44.4(2) 


-12.4(2) 


-2.88(2) 


0.176(3) 


1.29(1) 


4 He 


-105.4(4) 


-30.9(5) 


-16.3(1) 


0.63(1) 


7.26(7) 


6 He 


-127.(1) 


-44.(2) 


-20.3(4) 


-0.91(6) 


9.6(3) 


6 Li 


-150.(1) 


-38.(2) 


-19.8(4) 


-0.44(5) 


9.1(2) 


7 Li 


-178.(2) 


-54.(3) 


-25.6(6) 


-1.1(1) 


12.3(4) 


8 He 


-153.(1) 


-66.(3) 


-25.6(6) 


-4.0(2) 


13.3(4) 


8 Li 


-211.(1) 


-67.(2) 


-34.2(5) 


-3.8(1) 


17.4(3) 


8 Be 


-234.(2) 


-69.(3) 


-38.5(9) 


-0.9(2) 


18.3(6) 


7 n 


-10.11(9) 


-49.(1) 


-0.07(5) 


-5.4(3) 


1.9(1) 


8 n 


-12.0(1) 


—61.(1) 


0.31(9) 


-5.9(4) 


2.6(2) 



TABLE X. Isovector and isotensor energies (a^' T in keV) for isomultiplets. 



A 


T 


n 


AV18 


AV18/UIX 


AV18/IL1 


AV18/IL2 


AV18/IL3 


AV18/IL4 


AV18/IL5 


Expt. 


3 


l 

2 


1 


732(1) 


762(1) 


757(2) 


757(1) 


752(1) 


760(2) 


759(2) 


764 


3 


1 

2 


1 


723(2) 


753(1) 


748(2) 


747(1) 


746(1) 


751(2) 


750(2) 


764 


3 


1 

2 


1 


731(7) 


753(8) 


754(7) 


763(7) 


748(8) 


754(7) 


753(7) 


764 


6 


1 


1 


1068(4) 


1102(9) 


1149(6) 


1172(6) 


1141(5) 


1147(6) 


1179(6) 


1173 


7 


1 

2 


1 


1586(7) 


1565(7) 


1613(8) 


1588(7) 


1554(9) 


1609(8) 


1610(8) 


1644 


8 


2 


1 


1457(7) 


1488(5) 


1581(7) 


1622(8) 


1631(9) 


1670(9) 


1637(6) 


1659 


8 


1 


1 


1636(9) 


1672(7) 


1762(8) 


1810(6) 


1759(8) 


1837(7) 


1774(9) 


1770 


6 


1 


2 




251(18) 


275(13) 


293(13) 


266(16) 


287(18) 


280(17) 


223 


8 


2 


2 




158(4) 


170(4) 


180(5) 


180(5) 


188(5) 


175(3) 


153 


8 


1 


2 




135(7) 


135(12) 


143(8) 


99(12) 


148(10) 


156(12) 


145 



TABLE XL Various contributions to the isovector and isotensor energies (in keV) computed 
with AV18/IL2. The definitions of the contributions are given in the text. 



A 


T 


n 


K CSB 


vci(pp) 


V 1,R 


yCSB _|_ yCD 


Total 


Expt. 


3 


l 

2 




14(0) 


649(1) 


29(0) 


64(0) 


757(1) 


764 


6 


1 




16(0) 


1091(5) 


18(0) 


47(1) 


1172(6) 


1173 


7 


1 
2 




22(0) 


1447(6) 


40(0) 


79(2) 


1588(7) 


1644 


8 


2 




18(0) 


1528(7) 


17(0) 


59(1) 


1622(8) 


1659 


8 


1 




23(0) 


1686(5) 


24(0) 


76(1) 


1810(6) 


1770 


6 


1 


2 




166(1) 


19(0) 


107(13) 


293(13) 


223 


8 


2 


2 




136(1) 


6(0) 


38(5) 


180(5) 


153 


8 


1 


2 




141(1) 


4(0) 


-3(8) 


143(8) 


145 
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TABLE XII. Experimental and computed energies (in MeV) of "wide" or experimentally un- 
known states. The experimental widths (in keV) of the states are also given. 





AV8' 


AV18 


UIX 


IL1 


IL2 


IL3 


IL4 


IL5 


Expt. 


r 


5 He(r) 
B He(± _ ) 


-23.85(4) 


-22.47(9) 


-26.9(1) 


-27.7(1) 


-27.7(1) 


-27.4(1) 


-27.5(1) 


-27.4(1) 


-27.52(2) 


650 


-23.17(3) 


-21.9(1) 


-25.8(1) 


-26.5(1) 


-26.4(1) 


-26.3(1) 


-26.1(1) 


-26.0(1) 


-26.32(20) 


5500 


6 He(l+) 


-21.58(4) 


-20.2(1) 


-24.4(1) 


-24.7(1) 


-24.5(1) 


-24.2(1) 


-24.1(2) 


-24.1(1) 






6 Li(2+) 


-24.12(4) 


-22.7(1) 


-27.2(1) 


-27.9(1) 


-27.9(1) 


-27.7(2) 


-27.9(1) 


-27.8(1) 


-27.68(2) 


1700 


7 He(i-) 


-22.01(10) 


-20.8(2) 


-24.3(2) 


-26.6(2) 


-26.5(2) 


-26.3(2) 


-26.1(2) 


-26.3(2) 






7 He(§") 
7 Li(|") 


-20.81(10) 


-19.2(2) 


-23.2(2) 


-24.7(3) 


-24.4(3) 


-25.0(2) 


-25.0(2) 


-25.0(2) 


-25.92(30) 


2200 


-27.52(5) 


-25.7(1) 


-31.3(1) 


-32.3(2) 


-32.2(2) 


-32.0(2) 


-32.1(2) 


-32.3(2) 


-32.56(5) 


875 


8 He(2+) 


-21.39(8) 


-19.6(2) 


-24.1(2) 


-26.8(3) 


-26.6(3) 


-26.2(3) 


-27.2(3) 


-26.6(3) 


-27.82(5) 


630 


8 He(l+) 


-21.2(1) 


-19.6(2) 


-22.7(2) 


-26.0(3) 


-25.8(3) 


-26.2(3) 


-25.8(3) 


-25.8(3) 






8 Li(0+) 


-33.5(1) 


-31.3(2) 


-36.1(2) 


-38.4(3) 


-38.4(3) 


-37.2(4) 


-37.8(4) 


-38.3(4) 






8 Be(2+) 


-45.6(1) 


-42.7(3) 


-51.5(2) 


-53.6(3) 


-53.5(3) 


-52.4(3) 


-53.1(3) 


-53.2(3) 


-53.46(3) 


1500 


8 Be(4+) 


-38.7(1) 


-36.2(2) 


-44.9(2) 


-45.5(3) 


-45.4(3) 


-45.0(3) 


-45.4(3) 


-45.9(3) 


-45.10(30) 


3500 


8 Be(3+) 


-31.2(2) 


-29.3(3) 


-34.9(3) 


-37.2(3) 


-37.1(3) 


-36.9(3) 


-38.0(4) 


-37.2(3) 


-37.26(3) 
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TABLE XIII. Computed and experimental spin-orbit splittings in MeV. 









L 


S 


AV8' 


AV18 


UIX 


IL1 


IL2 


IL3 


IL4 


IL5 


Expt 


5 He 


l- 

2 


3~ 
2 


1 


l 

2 


0.68(5) 


0.6(1) 


1.1(2) 


1.2(2) 


1.3(2) 


1.1(2) 


1.4(2) 


1.3(2) 


1.20 


6 Li 


2+ 


-3+ 


2 


1 


0.86(6) 


0.8(1) 


0.9(1) 


2.2(2) 


2.2(2) 


2.4(2) 


1.9(2) 




2.12 


7 Li 


i- 

2 


3~ 
2 


1 


1 

2 


0.39(9) 


0.5(2) 


0.3(2) 


0.2(3) 


0.6(3) 


0.6(3) 


0.4(3) 


0.3(3) 


0.47 


7 Li 


5~ 
2 


7~ 
2 


3 


1 

2 


0.89(8) 


0.7(2) 


0.8(2) 


2.2(3) 


2.2(3) 


2.0(3) 


2.4(3) 


2.0(3) 


2.05 


8 Li 


1 + 


-2+ 


1 


1 


0.3(2) 


0.2(4) 


0.6(2) 


1.3(4) 


1.7(4) 


1.1(5) 


1.1(4) 


1.6(4) 


0.98 


7 n 


3~ 
2 


i- 

2 


1 


1 

2 


1.53(6) 


1.65(7) 


1.5(1) 


2.8(3) 


2.8(3) 


3.6(4) 


2.4(4) 


2.3(3) 





TABLE XIV. RMS point proton and neutron radii in fm. 





AV8 




AV18/UIX 


AV18/IL2 


AV18/IL3 


AV18/IL4 


Expt 




V 


n 


P 


n 


P 


n 


P 


n 


P 


n 


P 




1.66(0) 


1.82(0) 


1.59(0) 


1.73(0) 


1.59(0) 


1.74(0) 


1.60(0) 


1.74(0) 


1.59(0) 


1.73(0) 


1.60 


3 He 


1.85(0) 


1.68(0) 


1.76(0) 


1.61(0) 


1.76(0) 


1.61(0) 


1.76(0) 


1.61(0) 


1.76(0) 


1.61(0) 


1.77 


4 He 


1.50(0) 


1.50(0) 


1.44(0) 


1.44(0) 


1.45(0) 


1.45(0) 


1.46(0) 


1.46(0) 


1.44(0) 


1.44(0) 


1.47 


6 He 


2.06(1) 


3.07(1) 


1.97(1) 


2.94(1) 


1.91(1) 


2.82(1) 


1.99(1) 


2.97(1) 


1.99(1) 


2.96(1) 




6 Li 


2.50(1) 


2.50(1) 


2.57(1) 


2.57(1) 


2.39(1) 


2.39(1) 


2.44(1) 


2.44(1) 


2.38(1) 


2.38(1) 


2.43 


7 Li 


2.29(1) 


2.47(1) 


2.33(1) 


2.52(1) 


2.25(1) 


2.44(1) 


2.32(1) 


2.52(1) 


2.26(1) 


2.44(1) 


2.27 


8 He 


1.93(1) 


3.22(2) 


1.98(0) 


3.17(1) 


1.88(1) 


2.96(1) 


1.86(1) 


2.92(1) 


1.82(1) 


2.88(1) 




8 Li 


2.31(1) 


2.73(1) 


2.19(1) 


2.65(1) 


2.09(1) 


2.45(1) 


2.11(1) 


2.51(1) 


2.07(1) 


2.43(1) 




8 Be 


2.42(1) 


2.42(1) 


2.48(0) 


2.48(0) 


2.44(1) 


2.44(1) 


2.48(1) 


2.48(1) 


2.39(1) 


2.39(1) 




''n 




3.08(1) 




3.09(1) 




2.92(0) 




2.85(1) 




2.85(1) 




8 n 




2.98(0) 




3.03(1) 




2.92(0) 




2.88(0) 




2.88(1) 
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TABLE XV. Isoscalar and isovector magnetic moments, calculated in impulse approximation, 
in nuclear magnetons. 







T 


AV8' 


UIX 


IL2 


IL3 


IL4 


Expt 












Isoscalar 








3 He- 


- 3 H 


2 


0.408(0) 


0.405(0) 


0.403(0) 


0.402(0) 


0.404(0) 


0.426 


6 Li 


n 
u 


n 89^1 "\ 


U.oZl ^1 J 


n si 7(~\ \ 


n 81 (\(\ \ 


n 81 q('\ \ 


n 899 


7 Be- 


- 7 Li 


1 

2 


0.904(6) 


0.90(1) 


0.894(1) 


0.895(1) 


0.898(1) 


0.929 


8 B- 


- 8 Li 


1 


1.31(1) 


1.307(2) 


1.276(1) 


1.287(1) 


1.295(1) 


1.345 












Isovector 








3 He- 


- 3 H 


1 
2 


-4.354(1) 


-4.340(1) 


-4.330(1) 


-4.316(1) 


-4.331(1) 


-5.107 


7 Be- 


- 7 Li 


1 
2 


-3.92(1) 


-3.96(1) 


-3.93(1) 


-3.84(1) 


-3.96(1) 


-4.654 


8 B- 


- 8 Li 


1 


0.39(1) 


0.40(2) 


0.369(9) 


0.34(1) 


0.377(9) 


-0.309 


TABLE XVI. Quadrupole moments, calculated in impulse approximation, in fm 2 . 






AV8' 


UIX 


IL2 


IL3 


IL4 


Expt 




6 Li 




-0.27(8) 


-0.1(2) 


-0.32(6) 


-0.35(6) 


0.27(5) 


-0.083 




7 Li 




-3.6(1) 


-4.5(1) 


-3.6(1) 


-3.4(1) 


-3.9(1) 


-4.06 




7 Be 




-6.4(1) 


-7.5(1) 


-6.1(1) 


-5.4(1) 


-6.6(1) 






8 Li 




3.5(3) 


3.0(1) 


3.2(1) 


3.2(1) 


3.4(1) 


3.19(7) 




8 B 




7.1(3) 


8.2(1) 


6.4(1) 


5.6(1) 


6.6(1) 


6.8(2) 
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TABLE XVII. Expectation values of 8v and the net change in the binding energy due to the 



boost correction (in MeV) for three of the Illinois Hamiltonians. 







e 

Of 




ov - 


i /t /"* t r \ 






AIM O /TT 1 

AV18/1L1 


A "\ T~l O /TT A 

AV18/1L4 


A ^ T1 O /TT r 

AV18/1L5 


A ^ J-\ O /TT 1 

AV18/1L1 


A"\ 7"1 O /TT /I 

AV18/1L4 


A ^ j-\ O /TT r 

AV18/1L5 


3tj7 l + \ 
H (2 ) 


0.406(4) 


0.394(4) 


0.395(4) 


C\ AAA / O \ 

0.000(8) 


—0.001(8) 


—0.012(8) 


Hed ) 


n on A ( A \ 

0.394(4) 


0.382(4) 


0.386(4) 


—0.003(8) 


0.001(8) 


—0.005(8) 


He(0 T ) 


O 1 o/o^ 
2.16{2) 






— U.1(J(4J 


— U. 1)8(4) 


— U.1U(4J 


He(lT J 


2.81(9) 


2.93(8) 


o oo/o\ 

2.83(8) 


0.1(1) 


0.1(1) 


0.1(1) 


u He(2^) 






2.78(8) 






A 1 /l \ 

0.1(1) 


6t 


2.81(8) 


2.96(9) 


3.04(8) 


0.2(1) 


0.2(1) 


0.2(1) 


7t ■ / 3 — \ 

7 Li(f ) 


3.9(1) 


A A/1 \ 

4.0(1) 


3.8(1) 


0.1(2) 


0.4(2) 


0.3(2) 


7 Li(i ) 






4.1(1) 






0.4(2) 


7 Li(D 






4.1(1) 






0.5(2) 


8 He(0+) 


4.2(1) 


4.5(2) 


4.5(1) 


0.5(2) 


0.0(2) 


0.2(2) 


8 Li(2+) 


5.4(1) 


5.9(2) 


5.3(2) 


0.3(2) 


0.3(2) 


0.3(2) 


8 Li(l+) 






5.2(2) 






0.1(2) 


8 Li(3+) 






5.7(2) 






0.2(3) 


8 Li(4+) 


5.0(1) 




5.3(2) 


0.7(2) 




0.6(2) 


8 Be(0+) 


5.6(2) 


5.6(2) 


5.8(3) 


0.1(3) 


0.3(2) 


0.3(4) 


8 Be(l+) 






5.3(2) 






0.3(3) 


8 n(0+) 




1.25(4) 


1.12(5) 




0.61(5) 


0.57(6) 
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FIGURES 



FIG. 1. Energies of ground or low-lying excited states of light nuclei computed with the AV18 
and AV18/UIX interactions, compared to experiment. The light shading shows the Monte Carlo 
statistical errors. The dashed lines indicate the thresholds against breakup for each model or 
experiment. 

FIG. 2. Three-body force Feynman diagrams. The first, a, is the Fujita-Miyazawa, b is two-pion 
S-wave, c and d are three-pion rings with one A in intermediate states. 

FIG. 3. Energies computed with the AV18/UIX, AV18/IL2, and AV18/IL4 Hamiltonians com- 
pared to experiment for narrow states. The light shading shows the Monte Carlo statistical errors. 
The dashed lines indicate the thresholds against breakup for each model or experiment. 
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Fig. 2 (Pieper, et al.) 




